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Abstract 

We use a superspin Hamiltonian defined on an infinite-dimensional Fock space with positive definite 
scalar product to study localization and derealization of noninteracting spinless quasiparticles in quasi- 
one-dimensional quantum wires perturbed by weak quenched disorder. Past works using this approach have 
considered a single chain. Here, we extend the formalism to treat a quasi-one-dimensional system: a quan- 
tum wire with an arbitrary number of channels coupled by random hopping amplitudes. The computations 
are carried out explicitly for the case of a chiral quasi-one-dimensional wire with broken time-reversal 
symmetry (chiral-unitary symmetry class). By treating the space direction along the chains as imaginary 
time, the effects of the disorder are encoded in the time evolution induced by a single site superspin (non- 
Hermitian) Hamiltonian. We obtain the density of states near the band center of an infinitely long quantum 
wire. Our results agree with those based on the Dorokhov-Mello-Pereyra-Kumar equation for the chiral- 
unitary symmetry class. 
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1. Introduction 

It has been 50 years since Anderson wrote his seminal paper that came to define the problem 
of Anderson localization, a tight-binding model (in three-dimensional space) for noninteracting 
quasiparticles subject to a static disorder (random uncorrected on-site energies) [jj]]. The insight 
of Anderson was to realize that disorder, however weak, can have dramatic consequences on 
the nature of the eigenfunctions. In three dimensions, disorder turns extended states of the clean 
system into exponentially localized states when the characteristic kinetic energy of the extended 
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states, measured with respect to edge of the band, is of the same order or smaller than the char- 
acteristic energy of the disorder. 

For fermionic quasiparticles, this ratio can be tuned by changing the number of occupied 
single-particle states, i.e., by changing the chemical potential. The value of this ratio at which 
single-particle eigenstates are neither exponentially localized nor extended, if it exists, is called 
a mobility edge. The mobility edge realizes a quantum critical point, i.e., there exists a diverging 
length scale as it is approached [^|,[||5j]. On the metallic side of this continuous phase transition, 
disorder can be treated perturbatively thereby defining the diffusive regime. On the insulating 
side of this continuous transition, it is the kinetic energy that can be treated perturbatively 
A solid analytical grasp of the quantum criticality at a mobility edge remains an outstanding 
problem in Anderson localization. 

The existence and characterization of a quantum critical point at a mobility edge for fermionic 
single-particle states was argued on the basis of a scaling hypothesis combined with perturba- 
tive calculations (see Refs. for reviews). The tool of choice to implement the scaling 

hypothesis in a system of infinite size has been to encode the effects of weak disorder on the 
fermionic quasiparticles into an effective interacting field theory called a nonlinear-sigma model 
(NLSM). The space spanned by the fields of the NLSM, the target manifold, is determined 
by the intrinsic symmetries obeyed by the random microscopic Hamiltonian [|7|^|J^, ^o|jrT| , p^| | . 
These intrinsic symmetries are the presence or absence of the symmetry under time-reversal, 
spin-rotation, or charge-conjugation, respectively By demanding that the time-evolution 
generated by the random microscopic Hamiltonian is unitary and that it preserves the fermionic 
anticommutation relations, the target manifold of the corresponding NLSM belongs to one of 



the 10 classical symmetric spaces [ 1 3 1 in either their compact, noncompact, or supersymmet- 



ric incarnations [gp. The notion of 10 symmetry classes has thus emerged as a corollary to the 
scaling hypothesis, here implemented by 10 families of NLSMs. The local and global prop- 
erties of a target manifold fixes the number of independent coupling constants that enters in 
a NLSM. From this point of view, a (nonexhaustive) classification of the mobility edges for 
fermionic quasiparticles in one, two, and three dimensions [|l4j] into universality classes can 
be deduced from the curvature and the homotopy group of the 10 classical symmetric spaces 



[y,|ljl7|,|lj|ljg0£^ This construction of a 

NLSM starting from a random microscopic Hamiltonian relies on the existence of a diffusive 
regime [pT|]. However, mobility edges in d = 1, 2, 3 dimensions are usually far away from the 
diffusive regime. The predictive power of NLSM beyond the existence of mobility edges is there- 
fore severely limited by ones ability to study NLSM at their strong coupling fixed points. 

The validity of the scaling hypothesis has been tested with the help of nonperturbative numer- 
ical calculations on random tight-binding or network models. The hypothesis of one-parameter 
scaling in quantum wires as well as the existence of mobility edges in d = 1, 2, 3 dimensions 
for different intrinsic symmetries of the microscopic models is consistent with diverse numeri- 
cal calculations in systems of finite-sizes (see Refs. [10,5^ for reviews). In this context, random 
network models play an important role [||j3|,^J^p|@||0,|70,|9|5^] . Consider for 
example the plateau transition in the integer quantum Hall effect (IQHE) (for a review see Ref. 
J5l|]). Chalker and Coddington (CC) have proposed a random network model to study the effects 
of disorder on the highly degenerate single-particle states in the first Landau level of a two- 
dimensional electronic gas in a strong magnetic field when the cyclotron length is much smaller 
than the characteristic length over which the random potential varies significantly [|38||. The CC 
network model simplifies greatly the study of the plateau transition in the first Landau level. For 
example, the position of the mobility edge that separates two integer Hall insulating states in the 
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CC network model follows from a duality argument. Properties of the mobility edge in the IQHE 
can thus be simulated numerically for large system sizes and large statistical ensembles for the 
disorder with the help of the CC network model. In addition, the CC network model has been 
used as a regularization of the corresponding NLSM. Namely, the CC model has been mapped to 
a certain one-dimensional antiferromagnetic (super)spin chain, with the space of states being the 
tensor product of alternating highest and lowest weight infinite-dimensional irreducible represen- 
tations (irreps) of a Lie (super)algebra, also loosely called "(super)spins", and the Hamiltonian 
(an infinitesimal version of the transfer matrix of the CC model) built out of generators of the 
same algebra. The (super)symmetry of the spin chain was the same as that of the target space of 
the NLSM [ ^p3]M , p5|]^ , p7|pq ] . This was an important conceptual step in the sense that it pro- 
vided the hope that the flow to strong coupling of the corresponding NLSM could be captured at 
the level of the (super)spin chain in a way similar as the flow from the 0(3)-NLSM to the SU (2) 1 
Wess-Zumino-Witten (WZW) conformal field theory is captured by the Bethe Ansatz solution to 
the quantum spin— 1/2 Heisenberg chain [ |59| ] . Although, this hope proved ephemeral in the con- 
text of the plateau transitions (for a review see Ref. [Q7J), it lead to remarkable exact results when 
the same approach was used in the case of the two-dimensional Bogoliubov-de-Gennes (BdG) 
symmetry class with spin-rotation invariance but with broken time-reversal symmetry, and the 
corresponding network model [^|^iH|^^|. 

Unfortunately, this success story is the exception rather than the rule at the present. For exam- 
ple, consider the cases of spinless fermionic tight-binding models with either random real-valued 
or random complex-valued hopping amplitudes, whereby the hopping amplitude is restricted to 
connecting pair of sites with each site belonging to a different sublattice of a two-dimensional 
bipartite lattice. Such tight-binding Hamiltonians have the property that their spectra of energy 
eigenvalues are symmetric under charge conjugation, i.e., they are left invariant under sign re- 
versal whatever the realization of the random hopping amplitudes. This is most easily seen by 
changing the sign of the single-particle eigenfunction for all sites belonging to one sublattice. As 
this property is also shared by the Dirac Hamiltonian in quantum chromodynamics for a given 
configurations of the gluonic gauge fields, it is called a chiral symmetry (chiral-orthogonal sym- 
metry when the Hamiltonian is real-valued and chiral-unitary symmetry otherwise) In one 
and two dimensions, the band center in the chiral-unitary or chiral-orthogonal symmetry classes 
realizes a mobility edge separating two insulating phases very much in the same way as is the 
case with the plateau transitions in the IQHE. 

The study of localization in the chiral-unitary symmetry class in one-dimensions was pio- 
neered by Dyson [|&|]. Dyson found a diverging density of states on approaching the band cen- 
ter for a random microscopic Hamiltonian describing a single-particle hopping with a random 
amplitude between the nearest-neighbor sites of an infinitely long open chain. The band cen- 
ter realizes a quantum critical point that separates two insulating phases. This quantum critical 
poin t has been investigated using various methods such as a mapping onto a random XY model 
[p^], solving a recurrence relation |j67|, 58], solving a stationary Fokker-Planck equation [|69|], 
the Berezinskii diagram technique [|70[71], solving a one-dimensional diffusion process with 
reflecting and absorbing boundary conditions at the left and right ends |72j], respectively, solv- 
ing the evolution equation of a transfer matrix [[73|], a mapping onto Witten's supersymmetric 
quantum mechanics [ |74]|75| ], a mapping to Liouville quantum mechanics [[76]], and, at last, a 
mapping onto a spin-Hubbard-like model [ [77] , [78| , |79] ] . (We note here that in this case the spin-like 
model involves two sites with a highest and lowest weight representations of a Lie superalge- 
bra paired in a supersymmetric fashion. In this setting one has to solve the problem of the two 
couples superspins, analogous to the "spin addition" in quantum mechanics, albeit more com- 
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plicated due to the infinite dimensionality of the superspins.) At the band center, the typical 
wavefunction is not exponentially localized [J68Q, while the logarithm of the conductance is not 
a selfaveraging random variable [[73]], as would otherwise occur on both sides of the mobility 
edge . A weaker divergence of the density of states on approaching the band center survives in 
the two-dimensional chiral-orthogonal and chiral-unitary symmetry classes [ p4pl| , p0[pl] ] . This 
divergence is captured by effective field theories for weak disorder that include two-dimensional 
NLSMs with compact, noncompact, or supersymmetric target manifolds, respectively. These 
two-dimensional NLSMs can be regularized by two-dimensional random network models [^9|. 
In turn, the latter can be mapped onto spin-Hubbard-like chains, which, however, have been in- 
tractable to this date. One of the reasons is that a straightforward generalization of the methods 
of Refs. [HH,||||,|6H|@||@|80] leads to spaces of states on the sites of the chain that 
are not irreducible representations of the corresponding Lie superalgebra. The structure of these 
spaces from the point of view of the supersymmetry of the underlying physical problem is a so 
far unresolved question. 

A NLSM approach to the model studied by Dyson is not possible as there is no diffusive regime 
in one dimension. This difficulty can be remedied by considering thick quantum wires, which are 
also often referred to as quasi-one-dimensional wires. A thick quantum wire with a large number 
of transverse channels N ^> 1 appears in the limit where both the bare conductance and the mean 
free path i scale with N. This limiting procedure justifies the derivation of a one-dimensional 
NLSM with any one of the 10 classical symmetric spaces as the target manifold. In particular, 
one can derive a one-dimensional NLSM with the target manifold corresponding to the chiral- 
unitary symmetry class. In the same way as a two-dimensional NLSM can be regularized by a 
spin chain, a one-dimensional NLSM can be regularized by a zero-dimensional spin-Hubbard- 
like model with one or two sites. Solving such a spin-Hubbard-like model is a reasonable first 
step towards understanding the critical properties of the two-dimensional models in the chiral- 
unitary symmetry class. 

The purpose of this paper is to calculate the density of states in the close vicinity of the band 
center of a quantum wire belonging to the chiral unitary symmetry class for an arbitrary number 
N of channels by using a representation of the problem in terms of an effective zero-dimensional 
"superspin" model. As we have pointed out earlier, in this case the space of states is not obtained 
from irreducible representations of a superalgebra. However, the Hamiltonian of the effective 
model can be written as a bilinear form in the generators of a superalgebra, and it is this con- 
struction that we refer to as "the superspin approach". The density of states has been obtained by 
the superspin approach only for A = 1 [ [77|j7q , [79| | . 

The density of states of quasi-one-dimensional wires with arbitrary A in the chiral symmetry 
classes has been obtained using an entirely different approach [82 83J based on the Dorokhov- 
Mello-Pereyra-Kumar (DMPK) equation. This is an equation for the joint probability distribution 
of the Lyapunov exponents that are related to the eigenvalues of the transfer matrix [ p4| | along the 
wire viewed as a multi-channel scatterer. There it was shown that the density of states depends 
sensitively on the parity of the number of channels A. When A is odd, the Dyson singularity 
is recovered. When N is even, the density of states is controlled by random matrix theory up 
to multiplicative logarithmic corrections. This paper is devoted to uncovering this parity effect 
within the superspin formalism applied to the chiral-unitary symmetry class, the simplest among 
all 10 symmetry classes from the point of view of the DMPK equation [ |82| , |8~5| ] . Whereas the 
approach based on the DMPK equation is geometric [ |8^ , p7| ], the method to be presented below 
is algebraic. However, the DMPK equation is limited to quasi-one-dimensional wires, while our 
hope is that the algebraic insights gained in this work can be of use to solving spin-Hubbard-like 
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models related to two-dimensional NLSMs in the chiral symmetry classes. 

The band center of a superconducting quantum wire with broken spin-rotation symmetry is 
another example of a mobility edge separating two insulating phases [ |83| , |88| ] . Remarkably, the 
density of states displays the same Dyson singularity as a chiral quantum wire with an odd num- 
ber of channels [ |83| , p8| l . This suggests a deep connection between the corresponding symmetric 
spaces that has been partially explored from a geometric point of view in Refs. [ |8^ Jffl| ] . Hereto 
we hope that the algebraic method used in this paper could shed some light on this connection. 

As our computation of the density of states for a quasi-one-dimensional wire in the chiral- 
unitary symmetry class is quite involved, we now summarize the main steps leading to the de- 
pendence 

{(ery 1 |kT 3 (£r)| , iVodd, 
(1.1) 
(er)|ln(er)|, iVeven, 

of the density of states v{e) at the positive energy e measured relative to the band center whereby 
< et <C 1 with the mean scattering time r of order N 2 £/v F (£ is the mean free path and v F the 
Fermi velocity at the band center without disorder). The divergence of the density of states when 
N is odd is nothing but the Dyson divergence for the single chain random hopping problem. For 
even N, there is a multiplicative logarithmic correction to the power law predicted by random 
matrix theory. 

In Sec. H we define a quasi-one-dimensional random hopping model on the square lattice 
of length L = Ma and width Na with a being the lattice spacing and L 3> Na. We also 
provide a supersymmetric representation for the global density of states. In Sec. the density 
of states in the thermodynamic limit L — > oo, L S> Na, is recast as a certain expectation value 
in the ground state of an operator H acting on the direct product of a fermionic and a bosonic 
Fock spaces with a positive definite scalar product. The operator H is an infinitesimal version 
of the transfer matrix along the wire written in second quantized language, and, therefore, we 
call H the (effective) Hamiltonian, even though it is not Hermitian so that we need to distinguish 
its right and left eigenstates. Hamiltonian H is quartic in terms of the fermionic and bosonic 
operators spanning the Fock space on which it acts. Moreover, Hamiltonian H, although not 
Hermitian, possesses a certain degree of supersymmetry, which allows us to construct its ground 
sate and thus to compute the density of states of the underlying quasi-one-dimensional wire. 
Finally, Hamiltonian H depends parametrically on the energy scale e at which the density of 
states of the wire is to be evaluated and on the energy scale 1/r characterizing the strength of the 
static disorder in the wire. The construction of the ground state of H is done in Sec. [|, where it is 
expanded in terms of certain basis states and where we find a recursion relation for the expansion 
coefficients. We do not solve this recursion relation exactly, but we solve it and calculate the 
corresponding density of states approximately close to the band center er <C 1. The approximate 
solution is first presented in Sec. || for N = 1 and N = 2, and then for arbitrary even or odd N 
in Sec. ||. We summarize the important steps of our derivation and make a comparison with the 
DMPK derivation in Sec. 0. 

2. The quasi-one-dimensional random hopping model 

In this paper we reconsider the density of states for a single quantum particle in a thick quan- 
tum wire. We start with a quasi-one-dimensional tight-binding system of spinless fermions with 
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nearest-neighbor and next-nearest-neighbor random hopping matrix elements. On the square lat- 
tice the tight-binding Hamiltonian for N coupled random hopping chains takes the form 

N M 

= — y"n;ij C n,i c n+lJ + ^n;ij' c i,?: C 7i+2,j + n - c - ) ■ (2-1) 

i,J=l n— 1 

Here, the indices i and j label the iV chains while n labels the rows of sites in the direction along 
the wire. We impose the periodic boundary conditions along this direction, so that M + 1 = 1, 
etc. The spinless fermions are represented by operators that satisfy canonical anticommutation 
relations. The nearest-neighbor hopping strength t n .^ consists of a large uniform part tS^ and 
a small random piece St^^, whereas the next-nearest-hopping strength t^.^ is assumed to be 
purely random and small. That is, the chains are only weakly coupled by small random hopping 
amplitudes. Perpendicular to the chains there is no notion of distance, as every chain is coupled 
with equal strength to every other one. We shall consider the quasi-one-dimensional limit, where 
the length L of the disordered wire is much larger than its width Na, with a the lattice constant. 

For vanishing next-nearest-neighbor hopping, t' n .^ = 0, the Hamiltonian (2.1) reduces to the 
chiral random hopping model. This model is special because its energy eigenvalue spectrum ex- 
hibits a symmetry under charge conjugation, i.e., the Hamiltonian changes sign under the unitary 
transformation 

<i) - ("I)" (cn,i><i) • (2-2) 
Hence, the eigenvalue spectrum of the chiral random hopping model has a reflection symmetry 
around the band center. As we shall see below, because of this extra symmetry the states near 
zero energy have localization properties that are dramatically different from those of the states in 
other parts of the spectrum. 

Under the assumption of weak disorder, 

St n;ij: t 'n;ij < *> (2-3) 
it is legitimate to take the continuum limit in the direction along the chains. In the absence of 
disorder the ground state of each chain consists of a filled Fermi sea bounded by two distinct 
Fermi points. The low energy excitations, the plane waves with momenta near these two points, 
are commonly called left- and right-movers, respectively. Linearizing the spectrum about the two 
Fermi points yields a kinetic energy that is a first-order differential operator. We thus model the 
dynamics of the single quantum particle in a quasi-one-dimensional geometry by 
L 

•K c = J dx ^ (x) h c (x) tp (x) , 
o 

h c (x) = -ia 3 d x + a v (x) + a^x) + a 2 v 2 (x) + <7 3 v 3 (x), 

where we choose our units such that h and the Fermi velocity v F are one. The spinors ip(x) are 
2iV-component vectors, vAx) (jx = 0, 1, 2, 3) are N x N random Hermitian matrices, and a „ 
denote the three Pauli matrices and the 2x2 unit matrix acting on the left-right-mover degrees 
of freedom. Due to gauge invariance the potential v 3 (x) can be chosen to be zero in a system 
with open boundary conditions. In the continuum language the chiral symmetry is implemented 
by the interchange of left and right movers and is represented by 

a 1 h c (x)a 1 = —h c (x). (2.5) 

Chiral symmetry is thus only satisfied if the random matrices v and v 1 vanish, i.e., a chiral 
symmetric disorder potential is here necessarily off-diagonal with respect to the left-right-mover 



6 



degrees of freedom. In the presence of time -reversal symmetry, a 1 h* c a 1 = h c , the disorder po- 
tentials satisfy the additional symmetry constraints 

v 0&= v 0;ijt V t,ij = v l;ijt v 2& = v 2;ijt ( 2 -6) 

where i,j = l,---,N. The disorder potentials i; are assumed to be independent and Gaussian 
distributed with zero means and and with variances 

[W*K;«0OL = 2 9A X ~ x ') Mi + W-l)5 a 6 Jk ] , (2.7) 
respectively, where fi = 0, 1, 2 and the chain index i,j,k,l = l,---,N. The disorder strength is 
denoted by and /3 is the Dyson index with j3 = 1 (2) in the presence (absence) of time-reversal 
symmetry. 

We want to compute the Green function of the Hamiltonian (2A) in order to obtain its density 
of states at the energy e, which we assume to be positive, e > 0, without loss of generality. The 
two-point single-particle Green function is given by 



G(x, i, a\x' , i', a'; e + irf) = ( x,i,a 



1 



(2.8) 



h c -e-ir] 

where a, a' = L, R are left-right-mover indices, = 1, • ■ ■ , N are channel indices, the in- 
finitesimal regulator r\ is strictly positive, r\ > 0, and \x, i, a) denotes a position eigenstate. The 
global density of states v(e) for a given realization of the disorder is then obtained from 

L 



v{e) 



lim — 

v^o L 



(It 



1 N 

-Y 







E - 

i=l a=L,R 



Im G(x, i, a\x, i, a; e + irf). 



It is convenient to introduce the notation 



177 



-is 



(2.9) 



(2.10) 



The Green function (2.8 ) is an analytic function of icj in the upper part of the complex plane for 
any realization of the disorder. By analytical continuation of e > to the upper imaginary axis, 
\ui is strictly imaginary (ui becomes strictly positive). If so, one has to analytically continue ui to 
the lower imaginary axis in order to obtain the density of states (2.9). Translation invariance is 
restored in the disorder average of the Green function, i.e., one can define 



G{x 



x : iuj) = 



1 

N 



N 

E 



E 

1 a=L,fl 



G(x, i, a\x',i, a; iu>) 



(2.11) 



where the square brackets on the right-hand side denote disorder averaging. Correspondingly, the 
mean global density of states is 



If i / 

[v(e)] av = — / dx lim lim — Im G{x — x ; iu>). 



(2.12) 



The global density of states becomes self -aver aging in the thermodynamic limit L — * oo, so that 
the brackets on the left-hand side of Eq. (2.12) can be omitted. 

The single-particle Green functions (2.8) or (2.11) can be represented by functional path inte- 
grals as long as iw belongs to the upper part of the complex plane. Unless specified, we assume 
that u > 0. To perform the disorder average we make use of supersymmetry, so that 
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G(x-x';iu) = — 



where the action 



N 



i=l a=L,R 



(2.13a) 



5 = / dx L , 



is defined as the one-dimensions integral of the Lagrangian 

H=^(ih B + w)il> + ?(ih c + u)C 

and the partition function 

Z:= / Dip* Dip DC D£ e~ s . 



(2.13b) 

(2.13c) 
(2.13d) 

We have introduced, for any channel index i = 1 , • • • ,N and any left-right mover index a = 
L,R, the complex-valued integration variables £* a (x) and £ ia (x), here related by complex 
conjugation, together with the pairs of independent Grassmann fields ip* a (x) and ip i a (x). In 
the fu ncti onal integral formulation, the periodic boundary conditions imposed on the Hamilto- 
nian (2T) translate as follows: the bosonic fields obey periodic boundary conditions, but the 
fermionic fields obey antiperiodic boundary conditions in the x direction. The assumption to > 
guarantees the convergence of the Gaussian integral over the complex-valued integratio n vari- 
ables a (x) and £ s a (x). The action and the measure of integration in the field theory ( 2.13a ) 
contains supersymmetries that rotate boson and fermion fields into each other. These symmetries 
guarantee that the corresponding partition function is unity in every disorder realization: Z = 1, 
which simplifies the computation of averaged Green functions considerably. The ensemble aver- 
age over the disorder configurations in Eq. (2.13a) can be performed analytically by means of a 
cumulant expansion. With the Gaussian distribution (2.7 ), this gives the exact result 



x = 



. N 

-y 



i=l a=L,R 



Dip* Dip DC D£ iP ua (x)iPi a (x') e 



where the effective action 



La 

S e & = J dxL eS (x) 



(2.14a) 



(2.14b) 



(2.14c) 



is defined from the effective Lagrangian 

2 

+ E ^ Tr + (2//3 - 1)3 X] (*)■ 

(14=0 

Here, Tr (• • • ) refers to the trace over the channel index, we made use of the Hermiticity of the 
three random N x N matrix v 1 2 , and we also introduced the N x N Hermitian matrices O 1 2 
and their transposed 3oi2 through their matrix elements 



3 -ij = ^i,L^j,L + ^i,R^j,R + + C.hC.h- 

3l;ij = $i,B$j,L + ^th^ 3 ,R + £i,fl£/,Z + C /. £.,./,• • 
3 2 ;ij = 1 (<l>i,R^j,L ~ ^i,L^j,R + €i,Rtj,L ~ £*,l£j,r) 



(2.14d) 



with i,j = 1, ■ • • , N. 

We are going to demonstrate how we can use Feynman's transfer matrix method to trade the 
one-dimensional field theory encoded by Eqs. (2.14a), (2.14b), and (2.14c) with u> > for a 
zero-dimensional Schrodinger equation with a non-Hermitian superspin Hamiltonian. We will 
show that the non-Hermitian Hamiltonian possesses real eigenvalues and has manifestly unitary 
"imaginary-time evolution". We will also show how the density of states (2.12) can be extracted 
from an expectation value in the nondegenerate ground state annihilated by this non-Hermitian 
Hamiltonian. 



3. Mapping to a transfer Hamiltonian 



To perform the mapping onto a superspin Hamiltonian, we start from the classical Hamiltonian 



N 

,dcg,a ^?,deg,a ) 

£=1 deg=0,l a=L,R 

that can be read off from the Lagrangian ( [2.1 3c[ ) once AN coordinates 

%,de g ,a> « = !)••• deg = 0,1, a =1,2, 

have been chosen and by use of the AN momenta 

6L 



Pi 



5 % a ' 



i= IV, 



= 0,1, a = 1,2, 



(3.1a) 



(3.1b) 



(3.1c) 



whereby 6/5qi. a denotes the right derivative. We are thus interpreting x as the imaginary time. 
There is considerable freedom in the choice of the q's and hence of the p's. The most obvious 
choice is 



for which the momenta are 

Pi = {Pi,Q,Pi,l) , Pi,0 = (-Ci,L> ti,R) . Pi,l = {-^L^Ir) > 

with the channel index running over i = 1, • ■ ■ , N, With this choice, 



£ = p ■ q + Tr (lu3 q + iuo3 + + i«23 2 ) > 



where we have introduced the short-hand notation 

E E 



N 

,dcg,aQi,deg,a 

i=l deg=0,l a—L,R 



(3.2a) 
(3.2b) 

(3.2c) 
(3.2d) 



and, for any i, j = 1, • • • , TV, Eq. ( |2. 14d[ ) now reads 

3o-aj = ~Pi,l,Llj,l,L + Pi,l,Rlj,l,R — Pi,0,Llj,0,L + Pifi,Rljfl,R> 
= +Pi,l,B,Qj,l,L ~ Pi,l,LQj,l,R + Pi,0,Rljfi,L — Pi,0,L<ljfi,R> 



(3.2e) 



\Pi,\,R^3,\,L +Pt,l,i9j,l,JJ + Pi,0,RQj,0,L + Pi,0,L9j,0,R) 



Observe that the bilinear form 3q = {do-,ij) nas tne undesirable property that it is not positive 
definite. Hence, the quantum operator corresponding to Eq. (3.1) is seemingly ill defined as it 
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has an unbounded spectrum of eigenvalues from below upon quantization by which the classical 
p's and q's are replaced by operators that obey the canonical supercommutation relations 



LQi,deg,a 5 Pi' ,deg 



' Qi, deg, aPi' dc; 



g ,« 



(-) 



deg deg' 



Pi' ,deg' ,a' Qi,deg,a 



^i.i' ^ ol.cl' i 



(3.3) 



for any i, i' = 1, • • • , TV, deg, deg' = 0, 1, and a, a' = L, R. To cure this difficulty it is conve- 
nient to redefine the vacuum in the fermionic sector of the theory to be the Dirac sea. However, by 
the spin-statistics theorem, this medicine is inoperative in the bosonic sector unless one is will- 
ing to either give up Hermiticity [j78[] or the positive definite scalar product of the Fock space on 
which the quantum operator corresponding to Eq. (3.1) is defined \19\. In this paper, we choose 
the former route, so our quantum Hamiltonian will act on a space with a positive-definite scalar 
product, but will not be self-adjoint with respect to this product. In this setting, one will have to 
distinguish between the right and left eigenstates of the Hamiltonian. 
For example, consider the quantum theory defined by the assignments 



V'i.L-^/i. ll>i,L->-ji, C 

with the only nonvanishing supercommutators 



i.L 



{/-/;} {/.•/;} 



i,3 = 



+ 5t, 



,N. 



(3.4a) 



(3.4b) 



The corresponding Fock space $ R is obtained by action of the creation operators on the vacuum 
state, i.e., the state |0) that is annihilated by the right action of f it f i ,b i , b i with i = 1,2,--- ,N, 
i.e., 3r is the linear span of the product states 



n (//)"" (/7P ("ip (C'io>, 



i=i 

where 

n f i > n f. = M) n b i ' n b i 
We assume that the vacuum state is normalized to one, 

(0|0) = 1, 



0,1,2, 



(3.5a) 



(3.5b) 



(3.5c) 

whereby the dual vacuum state (0| is the state annihilated by the left action of fj, f},b\, b\ with 
i = 1, 2, • • • ,N. The identification (3.4) leads us to define the coherent states 



IV>il} R 
IV>l)r 



ft 



, l(iPr\ = 

^■ Jl |o>, M = 



L/R 



|0>, h {£ L \ 



Ai 



(3.6a) 



, N is implied, with the overlaps 



■Ai 



(3.6b) 



where summation over i = 1 , • 

We observe that |V>z,)r an d are not related by the "conventional" adjoint operation (• • • )T. 

Instead, the right coherent states |V , x)r can be mapped onto the left coherent states |?/> L ) L by a 
unitary transformation U, 



\^l)l = (MY =u\ip L ) 



R • 



U := cxp (mflfi 



(3.7) 
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Due to this fact, it is necessary to distinguish between the right Fock space 3r and its dual, the 
left Fock space tnat are related by the unitary transformation (3.7). That is, for a given basis 



|\E , m ) R of the right Fock space 3r, we define the corresponding basis |^ m ) L of the left Fock 
space 5 L by acting with U on the right basis, i.e., | v l' m ) L = U |^ m )rv The dual Fock space J£ 
is then the linear span of the basis L^ml = (I^jtJl) • Similarly, we will need to distinguish 
between left and right eigenstates of the transfer Hamiltonian H. 

Equipped with the definition ( |3.(i| ) for the coherent states, we a re now ready to de rive th e 
effective transfer Hamiltonian. It is given by replacing L in Eq. ( 3.1a ) with £ off from Eq. ( 2.14c ). 
After disorder averaging, the transfer Hamiltonian reads 



H = H 



H u = uTr J = w (N f + N f - + N b + N- b ) 



(3.8a) 



H v = E % Tr [J^ + (2/P - l)J„jJ] , 
with the number operators and disorder-induced operators defined by 



TV 



E*&> 



i=l 



N 

E^ 

i=l 



N 

E 

i=l 



A 7 * 



TV 

E/7/„ 

i=l 



(3.8b) 



and 





■■=+flfr 


-/<#+ 




Jl;ij 


■■=-fin- 








■■ i (/;./ 


1+fifj 


bP 3 + b,h) 



(3.8c) 



respectively. In as much as the right Fock space 3" R and its dual, the left Fock space that are 
associated with H are not related by the conventional adjoint operation, the transfer Hamiltonian 



H itself is not invariant under the conventional adjoint operation. On the one hand, Jq ^ = J ( 
On the other hand, because of the fermions, j| 4 • ^ J 1 ^, and y 7^ J 2 Consequently, 



if v 7^ if 



(3.9) 



Upon ma king th e transition to the transfer Hamiltonian if as encoded by Eq. (3.4) the partition 
function (2.13d) becomes 



Z=^ L ($ m |(-l) 



A F := JV) + A 7 /, 



(3.10) 



where Ap is the total fermion number operator, and the sum runs over some basis of the Fock 
space pair ($ L , S" R ), with |^ m ) L = ^l^m)^ f° r a ^ m - The factor (—1) n f is a consequence of 
the antiperiodic boundary conditions obeyed by the original Grassmann integration variables at 
the ends of the wire. 

Finally, we note that for a single channel, N = 1, alternative representations of Hamilto- 
nian (3.8a) have been studied in Refs. [ |77| | , f78| ], and [ |9l| ] as well as in connection with the 
random network problem introduced in Ref. [|79|]. We are now going to show that the eigenvalues 
of if for the chiral unitary symmetry class are manifestly real and that the matrix elements of 
cxp(— Lif ) between left and right eigenvectors of if define transition amplitudes of a unitary 
evolution operator, in spite of the fact that if is non-Hermitian. 
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3.1. Properties of the transfer Hamiltonian 



In this section and thereafter we shall limit ourselves to the chiral-unitary symmetry class 
defined by the conditions 

9o=0i = O, 13 = 2, (3.11) 



in Eq. (2.7). Hamiltonian (3.8a) thus reduces to 



H 2 , 



H„, = uTr J, 



H 2 = g 2 ^J 2 J 2 , 



(3.12) 



in the chiral-unitary symmetry class. The original supersymmetry of the action ( |2.14b| ) is realized 

ite 



by the invariance of Hamiltonian ( 3.8a[ ) under the interchanges 

fi ! fi i ~^~fi 5 fi 



i = l,---,N. 



(3.13) 



(3.14) 



In addition, for the chiral classes (g = g ± = 0), H 2 is unchanged by the interchanges 

(ft , fl+biM) "(fl,fi,-btK), i = l,---,N, 

whereas changes its sign under this transformation. This symmetry has its origin in the chiral 
symmetry. We also note that H 2 commutes with the eight generators of the Lie superalgebra 
gl{\\l) © <^(1|1), while commutes with the four generators of the diagonal sub-algebra 
C gl(l\l) © gl(l\l) (see Appendix^). 
The invariance of H under the supersymmetric transformation ( 3.13| ) has two consequences. 
First, it implies that an eigenvalue E l of H can only be nondegenerate if its right eigenvector 
belongs to a singlet, i.e., if the right eigenvector is annihilated by the supersymmetric transfor- 
mation ( 3.13 ). Second, any nonvanishing eigenvalue E l must be at least twofold degenerate with 
a pair of right eigenstates that are simultaneous eigenstates of the total fermion occupation num- 
ber N F that differ by an odd number of fermions. That is, all the eigenstates with nonvanishing 
eigenenergy E t can be grouped into supersymmetric multiplets of degenerate eigenstates. The 
only exception to this rule, is the nondegenerate ground state, a singlet state, with E l = 0. Since 
the transfer Hamiltonian H is non-Hermitian we need to distinguish left eigenstates from right 
eigenstates. In order to do so, we define right eigenstates, | ip l ( ) R G #r > with eigenvalues E t by 



E, 



(3.15a) 



and left eigenstates, L (<^ ; J € of H with eigenvalues E t by 

l(^J^ = l(^,J^ <=> H^\ v >u) L = E l \ V > u ) L , (3.15b) 

where the index 1 labels the different elements of the supersymmetric multiplet with eigenenergy 
E,. Since the right and left eigenstates are elements of the Fock space ^ R and respectively, 
it is possible to map \cp t t ) R onto \ip t t ) L by use of the unitary transformation U as defined in 
Eq. @, 

\<Pi,<,) h = U\<PtX lM=rWu\uI (3.16) 

Conversely, the unitary transformation U can be used to compute the adjoint of both J 2 and the 
transfer Hamiltonian H. Namely, we find that 

4 = U J 2 T U\ H ] = UHU\ (3.17) 

where (• • • ) T denotes the transpose operating on the N scattering channels. 
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As for the single channel problem of Ref. [|77|, we can now use the left and right eigen- 
states (3.15) to construct a normalized eigenbasis and the corresponding resolution of identity. 
Thereto, we argue that the left and right eigenstates of H in any given supersymmetric multi- 
plet with eigenenergy E l can be grouped into pairs (|</? ; t ) R / L , \ fi a^R/h) normalized by the 
condition 

l(Vi,i|¥Vm)r = 1 ' (3 - 18a) 
and such that the left and right eigenstates form a well defined biorthogonal system with 



fil,l'fi<T[L],L' 



(3.18b) 



for all energy eigenvalue indices I, I' and for all supersymmetric multiplet index i, if. Here, we 
have introduced the permutation function a[i] of the index i that parametrizes the pairing of the 
eigenstates. The permutation function is involutive, i.e., it satisfies a 2 = 1. The resolution of 
identity is then given by 

1. (3.19) 



E 



These properties of the eigenstates can now be used to show that the transfer Hamiltonian 
H with its associated Fock space pair (S^jSr) possesses a real eigenvalue spectrum, and that 
the operator exp(—LH) defines a unitary time evolution with respect to the Fock space pair 
3r)- First, for any given eigenenergy E l with left eigenvector \ [ <p l ^ and right eigenvector 
k, CTW ) R we find 



where we have used the mapping between left and right eigenstates (3.16), the Hermitian adjoint 
of H (3.17), and the normalization of the states ( 3.18| ). Second, the transition amplitude from the 
right energy eigenstate | (p l ( ) R to the left energy eigenstate | , ) L is defined by 



W, 



u,l' 



(3.21) 



after the analytic continuation L — > it. Using Eqs. ( |3.15[ ), Eqs. ( |3.18[ ), and the involutive property 
of a[i], one verifies that, for all t, 



El I w = Er<^ 

Z',t' l',i> 



■,+itH 



s +itE„ 



,-itE 



>6„,5 



V .1 <r[t'],t 



(3.22) 



I'd 



Hence, the expression ( 3.21 ) constitutes a well-defined transition amplitude, since it satisfies the 
condition of probability conservation assoc iated to a unitary time evolution. 

Summarizing, we have found that Eq. ( 3.17| ) implements the adjoint operation for the non- 
Hermitian Hamiltonian H. The operator U maps the right eigenstates from 3r into the dual 
space 5 L (i.e., the space of left eigenstates) in such a way that the transfer Hamiltonian becomes 
Hermitian within the Fock space pair (3l,3p{). That is, in order to reinstate unitarity of the 
evolution operator, exp(— itH), it is necessary to include the action of U in the scalar product. 
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3.2. Quantum representation of the density of states 



We are going to give a quantum representation of the density of states ( 2.12 ) in the long 
wire limit L — > oo and show how the original supersymmetry and a well-defined density of 



states (2.12) imply that the spectrum of H is positive definite. 



We start by re expressing the partition function ( 2.13d ) [see also Eq. ( 3.1C )] in terms of the 
eigenbasis (3.15) 



Z 



E 



-LE, 



1. 



(3.23) 



By the construction of Section |[ the supersymmetry implies that the partition function is exactly 
one whatever the length L of the wire. In the transfer Hamiltonian representation of the partition 
function (3.23), this is born out by the fact that each supersymmetric eigenmultiplet of H contains 
equal numbers of fermionic (N F odd) and bosonic (N F even) eigenstates, which thereby cancel 
in the supertrace due to the factor (— 1)^. Hence, the sum in Eq. (3.23) must then reduce to 
the expectation value in the zero-energy eigenenergy sector E Q = of the pair This 
sum then gives unity if and only if the eigenspace E = in # R ) is of dimension one. For 
the partition function to be independent of the wire length L, the spectrum of H must contain a 
nondegenerate eigenenergy E Q = 0. 



In the quantum representation, the density of states ( 2.12 ) is expressed by 



u{e) = lim. lim tt^Rc £ L ( Vu \ (-1)^ (B - B) 



— LH 



<pi,X 



lim lim TT^Re J^^M ("^^ " Q) ^ l^> H ' 



(3.24a) 



where either 



B 



N N 



(3.24b) 



JV 



E 6 ^ 



N l 



A'. 



(3.24c) 



As the system becomes infinite in length, L — > oo, any exponential term in Eq. (3.24), vanishes 



for all positive energy eigenstates. In order for the density of states to be well defined in the 
long wire limit, the eigenvalues of H need to be positive. If so, the density of states in the limit 
L — > oo is dominated by the ground state expectation value with E Q = 



v{e)= lim TT^Re^l (B-B) |p ) H 



L Re L (v? | (Q - Q) \<p ) R , 



lim 7r 

UJ — > — IS 



(3.25) 



where \(p ) h and \<P g) n d enote the nondegenerate left and right ground state wave functions, 
respectively. In Eq. (3.25) we have dropped the factor (—1)^ since the left and right ground 
state wave functions contain an even number of fermions, as we will explicitly verify in Section || 
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4. Construction of the ground state sector 



In order to compute the density of states in the long wire limit ( 3.25 ), we need to make an 
appropriate Ansatz for the solutions | ip ) R and | ip ) L of the right and left Schrodinger equations 

H Mr = (H u + H 2 ) |p ) R = 0, fft Mh = ( H l + h£) \<p ) L = 0, (4.1) 



respectively. The left and right ground state wave functions are related by |<^ ) L = U |v?o)r- 
As dictated by supersymmetry, the ground state wave function is nondegenerate and constitutes 
a supersymmetric singlet, i.e., it is annihilated by the supersymmetric transformation (3.13). 
Therefore, the ground state sector can be expanded in terms of states that transform as singlets 
under the supersymmetry (3.13). One such state is the vacuum |0). Other states that transform as 
singlets under the symmetry ( 3 . 1 3 ) can be generated by repeated action of H on the vacuum 1 0) . 
All the states that are generated by repeated action of H on the vacuum |0) form a vector space, 
which we call the right ground state sector R . The dual ground state sector 0£ is generated by 
repeated action of H< on the vacuum (0|. Alternatively, (5 L can be obtained by the application 
of the unitary transformation U defined in Eq. (3.7) on the right ground state sector © R . In turn, 
©l follows from © L with the help of the adjoint operation. The ground state sectors are subsets 
of the corresponding Fock spaces, i.e., <5 L C ^ L and 25 R C $ R . We shall thus expand the right 
(left) ground state \(p ) n (|<£>o)l) °f H ina basis of the right (left) ground state sector ©r (<3l)- 
Thereto we construct in what follows a right basis 



and a left basis 



{|m>W} 
{W(m|} 

^ J m.n 



u 



m, n = 1, 2, • ■ • , 



m, n = 1, 2, • ■ ■ , 



(4.2a) 



(4.2b) 



of the ground state sectors © R and (5 L , respectively. The left orthogonal set |to)l an d tne right 
orthogonal set |m) R "' ) are here generated with the help of the auxiliary orthogonal states 



|2m + 2,2n;|) := 



{A 



f n (D_) 2m+2 -2nS_S + (A_ 



2n-l , s2m+l 



10} 



|2m + l,2n+l;T) := [(A_) 2n+1 (D_) 2m+1 + (2m+l)S_S + {A_) 2n (D_) 2m \ |0) 
|2m,2n + 2;t> := [(A_) 2n+2 {D_f m + 2mS_S + (A_) 2n+1 (D_) 2m " r 
|2m+l,2n+l;4) := [(A„) 2 " +1 (£>_) 



. 2m+l 



(2n+l)S_S + (A_) 2n (D_) 2m \ |0), 

(4.3a) 



defined by taking arbitrary integer powers of the raising operators 

N N N 



N 



A_:=J2flfl D_:=J2blbl S_:=£/t&L ^ + :=^ & t /a t ; (4 . 3b) 



a=l 



acting on |0). The basis sets |w)l™' > and \m}^ l> are countably infinite due to the presence of 



>(n) 



boson creation operators. The basis set Im}^ is given by 
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N even 




(a) 

I 2m, N;T> 

1 2m+l, N-l;i) 
\2m+l, N-uT) 



Nodd 




(b) 

l2m+l,JV;t> 

1 2m, N-hl) 
\2m,N-l;1) 




\2m+l, 
\2m+l, l;T> 

1 2m, 0; i> 




2m+l, 
2m+l, l;T> 

2m, 0;i> 



Fig. 1. (Color online) These weight diagrams depict states (4.3; ) in the (N b , Nj,) plane. The left panel (a) displays the 
case of an even channel number N, whereas the right panel (b) shows the case of an odd channel number. 



in 



,(2n+l 
'R 

v(2n+2) 
' R 



where m = 1,2,3 

I \(2n+l) 

\m)l := 

I \ (2n+2) 

\m)l 
where m = 1,2, 3, 



) _ 1 (iV-2n-l) 



V2{N 


f 2m - 


-1)! 


1 (N 


- 2n- 


2)! 


V2(N- 


f 2m - 


-2)! 


, ■ • ■ . The basis 




1 (N 


-2n- 


1)! 


V2(N- 


f 2m- 


1)! 


1 (N 


-2n- 


2)! 



+ |2m- l,2n- 



^(\2m,2n;l) + \2m-l,2n+l;V 

( |2m - 2, 2ri + 2; |) + |2m - 1, 2n + 1; J.}) 



(4.4a) 



s/2 (iV + 2m-2)! 



is given by 

|2m,2n;|>- 
|2m - 2,2n 



|2m 



2;T) 



l,2n- 
- |2m 



i;T) 

l,2n 



L;I)) 



(4.4b) 



The states (4.3a) are eigenstates of the operators Nj, Np N b , and iVg defined in Eq. (3.8b), 
whereby the eigenvalues of 7Vy — TV^ and JVj, — JVg are vanishing for all the states in (4.3a). It is 
therefore convenient to depict these states in the (N b , A r ^)-plane (see Fig. |l|). As a consequence 
of Pauli's principle, the fermion number is restricted to a finite range. That is, any power 
larger than one of any of the fermionic operators S_ and S + vanishes, and for a given channel 
number TV we have 



N+l 



= o, 



and 



(A_) N D_ 



NS-S+ (A-) 



N-l 



(4.5) 



(4.6) 



Condition ( 4.5 ) determines the range over which the index n runs in Eqs. (4.3) and (4.4). We note 
that there is, as we shall see, a fundamental difference between even and odd channel numbers. 
For even the states with the highest fermion number N* are labeled by even boson numbers 
N b , whereas for N odd they are labeled by odd boson numbers N b . The right and left bases (4.4) 



are related by the unitary transformation U of Eq. (3.7), through 

U\m)g>, m = l,2,-.., 7i = l,2 r 



v(n) _ 



,N, 



(4.7) 



as it should be. The norms in the right sets (4.4a) are 
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(2n+l), 
R > 



v(2n+l) _ 
/R 



= (2m + 2n)(2m- l)!(2n)! 



(aT,+ H ) (m|m)g n+2) = (2m + 2n)(2m-2)!(2n + l)! 



(JV-2n-l)! 
(N + 2m- 1)1 
(JV-2n-2)! 
(7V + 2to-2)! 



=■ N« 



^(2) 



(4.8) 



with m = 1, 2, 3, • • • . It follows from ( |4.7| ), that the norm for a given state in the left set is equa l 
to th e nor m of the corresponding state in the right set. Moreover, the left and right sets ( A Ac ) 
and ( 4.4b ) are biorthogonal in the sense that 



{n l(m\m') { : ] 



to, to' — 1, 2, • ■ ■ 



1,2,"- 



(4.9) 



as follows from Eq. (4.8) and the orthogonality of the states Eq. ( 4.3a ). 

Equipped with a basis for the right and left ground state sectors, we are now in a position to 
expand the right ground state 



k>r = (r<^o|) =« (0) io) + EE 



n—1 m—1 



in terms of the basis Eq. ( |4.2a| ) and the left ground state 



AT oo 



ko>L = (i>oi) =« (0) io> + EE a » 



An) 



(4.10a) 



(4.10b) 



in terms of the basis Eq. (4.2b). The coefficients of these two expansions are identical and are de- 
termined by solving the Schrodinger equations (4.1). The overlap between the expansions (4.1C) 



for the right and left ground states follows form Eq. ( |4.9| ), 

l(</>oK)r = | a(0) | 2 - 
Consequently, in order to normalize the ground state wave function, we set 

a<°> = 1 



(4.11a) 



(4.11b) 



from now on. 



4. 1 . Ground state Schrodinger equation 



Hamiltonian ( 3.12 ) depends on two energy scales: the (imaginary) energy u at which the 
density of state s in the thermodynamic limit is to be evaluated and the chiral disorder strength g 2 
defined in Eq. (2.7). Since the density of states in the thermodynamic li mit is controlled solely by 
the right |</5 )r an d ^ l \ fo)h eigenstates annihilated by Hamiltonian ( 3.12 ), it must be a scaling 
function of the dimensionless variable 



U!, 



(4.12) 



provided we set g 2 to be unity, as we shall do from now on unless stated otherwise. The right 



(left) ground state |</>o)r (IVo)l) sat isfies the eigenvalue equation (4.1). In the basis (4.2), the 



eigenvalue problem (4. 1 ) yields the recursion relation for the coefficients a 
by 



(n) 



of Eq. (4.10) given 
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2(2m + 2n) uj a T (2n+1) = - 2(2m + 2n)(2m -l-2n + N)a^ n+1) 



(2n+l) 
m+1 



+ 2m(2m + l)a, 
+ (TV + 2?7i - 1){N + 2m - 2)a^"t 1} 
- (iV- 2n){N- 2n+ l)a^ n_1) 
-(2n + l)(2n + 2)a( f 2 l "+ 3 ), 
when n = 0, 1, • • • , [(N - l)/2] and 

2(2m + 277)wa^ n+2) = - 2(2m + 277)(2m - 3 - 2n 
+ 2m(2m-l)a%l + 1 2) 
+ (N + 2m-3)(N + 2m 



N)a 



[2n+2) 



(2n+2) 



m— 1 



(4.13a) 



2), 

-(N- 2n)(N - 2n - l)a<£ n) 
-(2n + 2)(2n + 3)aP n+4 >, 

when 77, = 0, 1, • ■ ■ , [N/2] — 1 and where m = 2, 3, 4, ■ ■ • . Here, it is understood that a 
when 77, < or tt, > N. The initial conditions are 

4w<4 X) = +V27V(yV + 1) - 4(7V + l)a[ 1] + 6a^ - 2aJ 3) , 

4waf } = -V2N(N - 1) - 4(iV - ljof 5 



(4.13b) 



2<4 2) -6<4 4) . 



(n) 



(4.13c) 



Equations ( 4.13 ) are a generalization to the case of an arbitrary number N of channels of the 
recursion relations obtained by Refs. [ |77|]78| , ^T| ] when N = 1. The rest of this paper is devoted 
to solving them and evaluating the density of states using Eq. (|3.25|). 



5. Evaluation of the density of states when N — 1 and N = 2 



Before we solve the recursion relations (4.13) for arbitrary N, we first specialize to the case 
N = 1 and N — 2 for pedagogical reasons. Only the case N = 1 using the supersymmetric 
quantum representation (see Refs. [ |77[[78| , ^T| |) was treated in the literature. The case N = 2 was 
solved with the help of the DMPK equation in Refs. [p2|] and | 



5.1. Density of states for a single channel N = 1 



For a single channel, N = 1, the recursion relations Eq. (|4.13[) for the coefficients dm simplify 



to 



2u>a{ 1} = +V2 - Aa\ L > + 3a 2 



(i) 



2ua 



(i) 

m+1 



-2(2m + 2)ai\ } +1 + (2m + 3)a£L + (2m + l)a«, 



(5.1a) 
(5.1b) 



with 777 = 1, 2, 3, • • • . This recursion relation is identical to the one obtained by Gogolin et 
al. @ by means of the Berezinskii diagram technique. Hence, one can identify the coefficients 



dm as the so-called right-hand side Berezinskii block. The recursion relation (5.1b) can also be 
written in terms of the finite differences 



AaL 1 ' 



a 



(i) 

m+1 



a 



(i) 



A 2 aL 1} 



a 



(i) 

rn+2 



2a 



(i) 

m+1 



a 



(i) 



(5.2a) 



That is, Eqs. (5.1b) and (5.2a) yield 
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2wo^ = +(2m + 3)A 2 a« + (2 - 2w) Aa&>, (5.2b) 
where m — 1, 2, 3 • ■ ■ . 

The strategy to analyze equations (5T) in the limit of asymptotically small frequency, w <C 1, 



is to derive two approximate solutions that are valid in two distinct asymptotic limits. For large 



m 3> 1, the variable m can be taken as continuous. This reduces the recursion relation ( 5. It ) 
to a Bessel differential equation. For small moj <C 1, we can neglect the u) term in Eq. ( 5.1b ). 
With this approximation, it becomes possible to solve the resulting recursion relation exactly, for 
example, by means of the generating function technique. 

For oj asymptotically small (the limit we are interested in), the overlapping region 1 -C m <C 
lu^ 1 becomes arbitrarily large. We can then derive matching conditions between the two ap- 
proximate solutions to fix the coefficients of the Bessel equation and the initial values of the 
zero-frequency recursion relation. Hence, in the limit of asymptotically small frequency, this 
procedure gives a complete solution for the ground state of H. This solution can be used to com- 
pute the expectation value of (B — B) in the ground state, which yields the asymptotic form of 



the density of states (3.25) as e — ► 



We first consider the limit m 3> 1. In this limit, we can neglect terms of order one compared 



to m in Eq. (5.2b) and replace finite differences by derivatives with respect to m, i.e., the discrete 



index m is now treated as a continuous one with a™ — > (m). This gives 

W 1 ' = 2 m-— + —— . (5.3) 
\ dm z dm J 

Here, we have assumed that a^ 1 varies slowly so that da^ / dm can be neglected relative to a^ 1 . 
After the substitution 

x 2 := Acorn, (5.4a) 

whereby 

J- » 1, (5.4b) 

this differential equation can be transformed into a modified Bessel equation and we obtain the 
solutions 

a«(z) = <Qk (x) + c { %(x), (5.5) 
where I (x) is the modified Bessel function of the first kind and K (x) the modified Bessel 
function of the second kind (see chapter 9.6. in |^|). We shall demand that a™ decays to zero 
as 77i — * oo, i.e., we must set 

£1 = (5.6) 



in Eq. (|5.5|). The remaining coefficient c^l of the modified Bessel function K is determined by 

matching the solution (5.5) to the one for m <C to" 1 . The coefficient c^l of the modified Bessel 
function K will thereby acquire an to dependence. In order to carry out this matching procedure, 
we will need the small x behavior of the continuum solution 



a W (x) ~ -c$ [ln(z/2) + 7 ] = — |i pn(mw) + 2 7 ] , (5.7) 

where 7 is the Euler's constant. 

Second, we consider the limit 777 <C in which case we can safely drop the terms contain- 



ing u> in Eqs. (5.1). Equation (5.2b) can then be written as 
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= (2m + 3)Ao^ 1 -(2m+l)AaW 
= A((2m + l)Aa£ 



with m = 1, 2, 3, • • • , together with the initial condition 



3a 



in 



4a 



(l) 



V2. 



(5.8a) 



(5.8b) 



Integrating Eq. (5.8a) and determining the integration constant by use of Eq. (5.8b), we find the 
first order difference equation 



(2m+l)(a« 1 -a«)=4 1) -^ 



(5.9) 



with m = 1, 2, 3, • • • , whose solution is given in terms of the Digamma function (for refer- 
ences on the Digamma function "J, see chapters 6.3 or 8.36 in Refs. [^2|] or [^3||, respectively) 



a® =^2 + \ (a« - V2) (*(m + 1/2) - *(l/2)) , 
with the limiting form 



a W ~ V2+ 1 (a« - V2) (lnm- *(l/2)) 



(5.10) 



(5.11) 



for large meN. 

Having solved the recursion relation in the two limits 1 <C m and m <C oj^ 1 , we now match 
solutions ( |5.7| ) and ( 5.11 ) in the overlapping regime 1 <C m -C ijJ -1 to fix the coefficient 

of the modified Bessel function K and the initial value of the u> = recursion relation. 
Equating Eq. (5.7) with Eq. (5.11) yields the two conditions 



-o,l 



V2, 



-eg Inc. - 2c« 7 = 2V2 - (a« - V2) *(l/2), 



(5.12) 



for the two unknowns c^ 1 ] and . Solving for and ajj 1 , we get 



V2 



(ln(C» 



-2\/2 
In (Cjw) 



, C 1 :=e«p(27 + *(l/2)). (5.13) 



We have thus derived two limiting solutions to the recursion relations (5.1 ) 



^(l 



^(m+l/2)-^(l/2) 
lnfCjw) 



1 < m, 



(5.14) 



Knowledge of the asymptotics (5.14) is sufficient to determine the leading behavior of the 
density of states v(e), Eq. (3.25), for asymptotically small energies e. The density of states in 
the long wire limit Eq. (3.24) is given by the expectation value of the operator (B — B) in the 
ground state 

v{e)= lim 7r _1 Re L (^ | (B-B) \ip ) R 



lim 

uj — > — ie 



I 2 \ 



(5.15a) 
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where we have used Eq. (1.8) and the identity 



(B-B)\ 



>(i) 



(5.15b) 



In view of Eq. (5.14), it is necessary to break the sum over m in Eq. (5.15c) into two parts 
separated by the integer m , with 1 <C m <C We choose m to be some fixed number, 
independent of us. With this choice and for sufficiently small u>, the sum on the right-hand side 



of Eq. (5.15a) is dominated by the part m > m , i.e., by the contribution from the first line in 



Eq. (5.14). (In the limit uj — > 0, the sum over m < m , i.e., the contribution from the second line 



in Eq. (5.14), is less divergent than the sum over m > m , and can therefore be neglected.) We 
then find the estimate 



M(B-B)\(p ) R 



ln 2 (C» 

4a;" 1 
ln 2 (C» 



dm Kn (2\/mLu) 



dx xKq (x) 



This gives the asymptotic behavior of the density of states 

1 

Re — ^ 

uj — > — is 



v{e) oc lim Re ■ 2 

win (Ciw) 



e \\tl e 



3 „| ' 



which has precisely the form of Dyson's singularity. 
5.2. Density of states for two channels N = 2 



(5.16) 



(5.17) 



For the case of two channels, N = 2, the recursion relations (4.13) for the coefficients a 



(2) 

and f reduce to 



,(!) 



with m = 1, 2, 3, 



2wa^ ) +1 = (2?n + 1 
and with the initial conditions 

2u}a i i ) = + V23 - 6a{ 1] 



2a (1) - 

zu m+l 


r u m+2 




2a (2) - 

zu m+l 


ha (2) - 





2a;a 



(2) 



y/2- 2a[ 2) + 



,( 2 ) 



(5.18a) 



(5.18b) 



Alternatively, we can rewrite the two decoupled recursion relations ( 5.18a ) in terms of the finite 
differences 



(0 

J m+1 



A 2 a« 



a W - 2a W 

u m+2 iu m+l 



1,2. 



With this, Eq. (5.18a) becomes 

2c a W=(2m + 3)A 2 a W-2 W A a (i> 

2wo« = (2m + l)A 2 a£> - 2wAa£> 



(5.18c) 



(5.18d) 



To analyze these equations, we follow the same strategy as in the previous section for the N = 1 
case and solve the Eqs. (5.18) first in the m ^> 1 limit and then in the m -C limit. 
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We begin with the m ^> 1 limit. Treating m as a continuous index, we replace finite differences 



by derivatives in Eq. (5. 18d ), i.e., we let a 



dm 1 



a (m) with i = 1,2. This gives 
i = 1,2. 



(5.19) 



Here, a^ 1 ' is assumed to be slowly varying so that we can neglect da^ /dm compared to 



After the substitution 
whereby 



4w 



4u>m, 



> 1, 



(5.20a) 



(5.20b) 



Eq. (5.19) reduces to a Bessel-type equation whose solution is given in terms of the modified 
Bessel functions K^x) and Ii(x), 



a«H*)=c^K 1 (x) + c^I 1 (x), 



i=l,2, 



(5.21) 



(i) 

with the coefficients j with i = 1 , 2 of the modified Bessel functions K 1 when j = 1 and I I 
when j = 2. We shall demand that a m with i = 1, 2 decay to zero as m — > oo, i.e., we must set 



= 0. 



(5.22) 



in Eq. ( |5.21[ ). The remaining coeffici ents with i = 1, 2 of the modified Bessel function K 1 
will be fixed by matching solutio ns (5.21 ) to the small mui solution. In order to do so we will 
need the small x behavior of Eq. ( 5.21 ) 

a«W~^«+cg^[ln(x 2 /4) + 2 7 -l] 



with i = l,2 and where 



C§° :=exp(2 7 -l). 



(5.23a) 



(5.23b) 



Next, we turn to the limit rruv <C 1, in which we can neglect the lo term in Eq. (5.18d). If so, 



Eq. ( 5. 18d) becomes 



together with the initial conditions 



2.(1) 



= (2m + 3)A 2 a, 
= (2m + l)A 2 a(, 2 l ), 



2<4 1} = +V2 

-V2- 



2a 



(2) 



„( 2 ) 



The solution of Eqs. (5.24a) is given by 

a$ =ma^ 1) - (m- 1)V2 , 
/2 



= m Oi + (m — 1 



(5.24a) 



(5.24b) 



(5.25) 



where we have fixed two initial values by use of Eq. ( |5.24b ). 

Having solved the two decoupled recursion relations (5.18) in the two limits m 3> 1 and 
m<w _1 ,we now match solutions (5.21 ) and (5.25) in the overlapping region in order to fix the 
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coefficients cj?\ of the modified Be ssel f unctio ns Kj and initial values with i = 1 , 2 of the 
to = recursion relation. Equating ( 5.23 ) with ( 5.25 ) and matching equal powers of m (whereby 
we neglect In m compared to In uf) gives the two equations 



W2 = I$ 



(5.26a) 



for the two unknowns c^l and and the two equations 



(2) 



v z — 2 C e,l : 

-e,i7T 



V5 = 



In G 



.(e), 



(5.26b) 



(2) (2) (1) (2) 

for the two unknowns c e { and a\ . Solving for a\ and a\ , respectively, we obtain 



,(!) 



(2) 



72 
n/2 



win (Cl e) w 



l + wlnfc^w) 



We have thus derived two limiting solutions to the recursion relations (5.18) 



(5.27) 



,(2) 



m \p2 lo In 



Km, 



(5.28) 



- 2V2mw ^ (2 v / mw) , 
— m 



V2, 



1 < m, 



In the long wire limit, the density of states is given by the expectation value of the operator 
(B - B) in the ground state [see Eq. fl3.24| )] 



v[e) = Urn. 7r _1 Re L (yj | (£? — S) |<^ ] 



= lira 7r 1 Re 



OO 

■E 

m— 1 



(a£>) (a^y 
2m + 1 2m - ] 



(5.29) 



Here, we have used Eq. (^J), the identity (B - B) |0) = 2 |0) and 



[B - B) \m 



[B - B) \m)£ =-\m 



A2) 
'B 



.(2) 



(5.30) 



We want to estimate the leading behavior of the density of states for two scattering channels, 
N = 2, and for asymptotically small energies e, i.e. uj <C 1. Thereto, we break the sum over m in 



Eq. (5.29) into two parts separated by the integer 1 <C m 



' 1 . With the choice tuq 



we find that the sum in Eq. (5.29) is dominated by the small m <C 1/lu solutions in Eq. (5.28) 
[the rest of the sum is easily shown to be 0(u>)], 
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1/u 



771 — 1 



1 



1 



2m + 1 2m - 1 



1 + mwln ( Cq^uj 



(5.31) 



This sum over m can be computed exactly by means of a telescopic expansion. We find 

4 



>ol(£-£)l¥>o>i 



2 + cj 
21n4-* 



2 - 2 7 



(w- 1 + l/2) -*(w~ 1 + 3/2) ]n(c£ c) w) 



(5.32) 



To compute the leading behavior of the density of states it is sufficient to retain only the lowest 
order in ui. This gives 



L ((^ | (B — B) |</3 )r ^ — w m2 w + higher order terms. 
We thus find the estimate 

v{e) oc lim it~ 1 Rc (— lo\yi 2 m) cxe|lne| 

uj— > — is 

for the density of states in the thermodynamic limit. 



(5.33) 



(5.34) 



6. Evaluation of the density of states when AT = 1, 2, 3, 



We are going to compute the density of states v{e) from Eq. (3.25) to leading order in the 
positive dimensionless energy e <C 1 by solving the recursion relation (4.13) for an arbitrary 



number N of channels. We begin with the case of N odd and proceed with the case of TV even. 
In doing so, we are going to reproduce the parity effect ( |1.1[ ) in the density of states that was 
obtained for the first time using the DMPK approach in Refs. [B2fl and [K3J. 



6.1. Density of states for an odd number N of channels 

The calculation of the density of states for an odd number N of channels fo llows along the 
lines of the N = 1 case with the caveat that for N > 1 the recursion relations ( 4.13 ) no longer 
decouple. This difference complicates the calculation considerably. To overcome this difficulty 

(n) 

we introduce a linear transformation of the coefficients a; n that approximately decouples the 



recursion relations (4.13) in the two regimes m> N and m <C ui 1 . That is, we introduce 



(N-l)/2 

e +i) = e [v,it«i 2n+i) , *=o,i,-, 

(JV-3)/2 

e +2) = E K ;J v, 2 ]!ai 2 " +2) , * = 0,l,- > 



N-l 

N - 3 
~2~' 



(6.1a) 



with m = 1, 2, 3, • • • and the transformation matrices 
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k (2n)!(7V-2n-l) 

L JVL o;AmJ„ —I 1 J „|fJV-l ^1 



3-P2 



1 1-iV 
2 ' 2 



(2n+l)!(AT-2n-2)! 
n!(^2-n)! 



3-P2 



-fc, fc + 2, -n 



3 3-JV 

2 ' 2 



;i 



(6.1b) 



(n) (k) 

It is also possible to express the a m s in terms of the 6 m s, 



Here, the symbol 3F2 denotes a generalized hypergeometric function (see for example Ref. J93|l). 

ess the ah 

(N-l)/2 

[ JVL o;iV : lJ °m ' ft — U, 1, , 2 ' 

(6.1c) 



e +i) = e 

fe=0 
fJV-3)/2 

E [ M o;]v,2 



(2n+2) 



fc=0 



& (2fc+2) 



77 = 0, 1, ■ • • 



2 ' 

AT-3 



(6.2a) 



With these definitions, we are going to rewrite the recursion relations (4.13) in a form that de- 
couples both in the limit uj = and in the large m limit. 

To see this, we first make use of identities (B.3a) and (B.3b) to combine Eqs. (4.13) and (5.1) 
into the recursion relations 

477iw6 1 ( r 2 l fe+1) = [-4m(2m — 1 + N) + (N — 1- 4k 2 )] b<£ k+1) 

+ 2m(2m + l)b% k + 1 1) + (2m - 2 + N)(2m + N- l)b { ^ ] 

(N-l)/2 

-** E "K/vJH 2n+1) 

with fc = 0, 1, • • • , (N - l)/2 and 

Amub^+V = { - 4m(2?7i - 3 + N) + [3N - 5 - 4(fc + l) 2 ] }6^ fe+2) 

+ 2777(2777 - l)6^+j 2) + (2777 - 3 + N)(2m -2 + N)b ( ^X 2) 

{N-3)/2 

-4a; E "K ; /v, 2 ]"^" +2) 



(6.2b) 



with k = 0, 1, ■ • ■ , (AT — 3)/2 and m = 2, 3, ■ • ■ . Here, the initial conditions are 

4ivb ( i k+1} = + V2N(N + 1) - (3JV + 5 + 4/« 2 )^ 2fc+1) + 66^ +1) 

(iV-l)/2 

-4c E n [ M o;N,i] k n a[ 2n+1) 

n=Q 

with jfe = 0, 1, ■ • ■ ,(iV-l)/2and 



(6.2c) 



Aujb 



(2fc+2) 



V2N(N - 1) - [N + 1 + 4(fc + l) 2 ]6f fe+2) + 2&!, 2fc+2) 

(JV-3)/2 

4. e -K^ 2 ]M 2n+2) 



(6.2d) 



n=0 



with fc = 0, 1, • ■ ■ , (AT — 3)/2. We aim at a solution of Eq. (6.2) for asymptotically small fre- 
quency cu. As in Sec. |5.1[ , our strategy is to solve Eq. ( |5.2] ) in the two limits N m and m <C 
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ui 1 . In the former limit, m ^> N, we treat m as a continuous variable. Then, Eq. (5.2) decou- 
ples, provided we assume that a~m decays rapidly for m ^> N. In the other limit, m <C lo^ 1 , we 
can neglect w in Eq. (|5.2|), which again decouples. In this way, it is possible to find approximate 
solutions in the two regions N <C m and m <C oj^ 1 that are uniquely fixed up to some mul- 
tiplicative coefficients and initial values, respectively. For asymptotically small w, the overlap 
between these two regions N <C m <C lo^ 1 is large. We can then match the two approximate 
solutions in the overlapping region. This gives a unique and approx imate solution for the ground 
state wave function, which in turn determines the density of states (3.25) in the long wire limit. 



6.1.1. Solution when m ^> N 



(n) 

First, we treat the limit m ^> N. If we assu me th at a,„ deca ys rapidly for large m, we can 
drop the last line on the right-hand side of Eqs. ( 6.2a ) and ( 6.2b[ ), respectively. Consequently, in 
terms of the finite differences 



AftW : = b . 



b (i) 



A 2 feW — 7>W 
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with i = 1, ■ • ■ , JV, the recursion relations Eqs. ( |6.2a ) and and ( |6.2b| ) read 

(m + l)4w (Afo( f f +1 > + C fc = - 2(2m + 2)(iV - 2)A6^ fe+1 ' 

+ (2m + 2) (2m + 3) A 2 b£ k+1 ' 
+ (N - 1){N - 2)b£ k+ V 

+ (N — 1 — 4fc 2 ) (A&^ fc+1) 
with A; = 0, 1, • • • ,(iV-l)/2and 
(m+l)4 W (A6^+ 2 '+6(f+ 2 



6« 



(6.3) 



(6.4a) 



6 (2fc+i; 



- 2(2m + 2)(N - 2)Afc< 2fc+2) 

+ (2m + 2)(2 7 n+l)A 2 fe( 2fc+2 ) 
+ (N ~ 2){N ~ 3)b^ k+2 ^ 

+ [37V - 5 - 4(fc + l) 2 ] (A6^ fc + 2 ) + C fe + 2 )) 



(6.4b) 



with k = 0, 1, • • • , (iV — 3)/2. In the limit m 3> TV, we can neglect terms of order iV compared 
to m and replace finite differences by derivatives. In place of Eq. (5.4) and if we assume that 
bin — > b (to) with i = 1,2, ■ ■ ■ , N is slowly varying, we get 



4mc^ 2fc+1 > 
with = 0, 1, 

4mwfe (2fe+2) 
with fc = 0, 1, 

whereby 



= -4m(N-2) 



db (2k+l 



dm 



4m 



d 2 & (2fc+l) 



dm 2 



[(N - l) 2 - 4fc 2 ] 6< 2fc+1 ) (6.5a) 



,(iV- l)/2and 



-4m(iV-2) 



-Am 



id 2 b (2k+2) 



dm dm 2 
, (N — 3)/2. By use of the substitution 

x 2 := 4wm, 



4^ 



+ [(iV-l) 2 -4(fc + l) 2 ]fe (2fc+2) (6.5b) 



(6.6a) 
(6.6b) 
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we find that the solutions to Eq. (6.5) are given by the linear combinations of 



with k = 0, 1, • • • , (N - l)/2 and 



JV-l 



K 2k {x) + c} o2 



(2fc+l) (X 



N-l 



(2/C+2) (X 



N-l 



K 2k+2( x ) + 4,2 



(2fc+2) fX 



N-l 



I, 



2fe+2 



(x) 



(6.7a) 



(6.7b) 



with k = 0, 1, • • • , (N — 3)/2, of modified Bessel functions if 2A . and 7 2/c . We shall demand that 
bm with i = 1, ■ • ■ ,N decay to zero as m — > oo, i.e., we must set 



c (l) - 

' 0.2 — u 



(6.8) 



with i = 1, • • • , N in Eq. (|677j). The remaining N coefficients c with i = 1, • • • , iV of the mod- 
ified Bessel functions i^ , i^ 2 , i^ 4 , • • • , K N _ 1 are fixed by matching solutions (6.7) to the solu- 
tions in the m <C uj^ 1 region. Thereto, we need to extract the terms that are of order x N+2k ~ 1 
and x N+2k ~ 3 from the expansion o f Eq. ( 5.7a ) and the terms that are of order x N+2k+1 and 
x N+2k ~ 1 from the expansion of Eq. (5.7b), when x <C 1. This gives 



& (1) (*)~cW(f) 



N-l 1 



- -2 7 - In — 



(1)/ iV^l 1 / (q) \ 

& {{moj) 2 - In ^6q 'tow J 



(6.9a) 



6( 2fc+1 )(x)^4 2 1 fe+1) 



„(2fc+l) 



(x/2) 



N+2k-l 



2(2k)\ 
(mcu)~ ' 



- 7 -In— +*(2jfe+r 



(x/2) 



Af+2fc-3 



2(2fc - 1)! 



l -+k 



2(2fc)! 



with ft = 1, 2, • • • , (N - l)/2, and 



„(2fc+2) 



(ar/2) 



JV+2A-+1 



2(2fc + 2)! 
(mw) 2 



In ( C { °\ 



—7 — In ■ 



(mw)" 



-+fe-i 



2(2ft- 1) 



(6.9b) 



*(2ft + 3) - 



(x/2) N + 2k - x 
2(2fc + 1)! 



In C 



(tow)" 



2(2fc + 2)! 
with ft = 0, 1, • • • , (N - 3)/2, and where 

C { ° } := cxp (7- #(2ft+l 



2(2ft+l)! 



(6.9c) 



(For references on the Digamma function "P, see chapters 6.3 or 8.36 in Refs. 
respectively.) 



(6.9d) 
or [61, 



6.1.2. Solution when 1 < m <C W 1 

Second, we treat the limit 1 < m <C w~ 



in which case we can neglect the uj terms in 



Eq. (6.2). In doing so, Eq. (6.2) becomes 

= [-4to(2to - 1 + N) + (N - 1 - 4k 2 )] b£ k +1) 
+ (2m - 2 + N)(2m - 1 + AO&m-l^ 



2m(2m + 1)6, 



(2fe+l) 
1 



(6.10a) 
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with/c = 0,l,--- ,(N- l)/2and 



(6.10b) 



= { - 4m(2m - 3 + iV) + [3N - 5 - 4(fc + l) 2 ] }6^ fc+2) + 2m(2m - 
+ (2m - 3 + N) (2m - 2 + iV)6^^ 2) 

with k = 0, 1, ■ • ■ , (AT — 3)/2, m = 2, 3, • ■ ■ , together with the initial conditions 

= + V2N{N + 1) - (5 + 3N + Ak 2 )b {2k+1) + 6b (2k+1) (6.10c) 
with A; = 0, 1, • • • , (AT - l)/2 and 

= - V2N(N - 1) - [AT + 1 + 4(fc + l) 2 ]b[ 2k+2) + 2b {2k+2) (6.10d) 



with k = 0, 1, • • • , (N — 3)/2. The solution to the recursion relation (6.10) can be expressed in 
terms of generalized hypergeometric functions 

6(1 , _ g(o,i) y 1 (m-l-i + ^l)! e(o , 2) (m-l + f)! 
m f-(2Z + l)(m-l-0! (m-i)! ' (6 ' Ua) 



u(2k+l) _ p(o,l) { m + k 2 



+ f) 



3^2 



+ g(o,2). 



(m - 1)! 
m + fc-l + f)! 



i — fc, 1 — k, 1 — m 
H-k-m-% 



;i 



"3-^2 



ft, 2 



;i 



with ft = 1, 2, • • • , (2V - l)/2, and 
(to + k — 



b (2k+2) = g 



(o,3) 



(m-1)! 



"3-^2 



— i — ft, — 1 — ft, 1 — m 



a(o,4) 



1 + f)! 



I)! 



1 I _ jU 

2 ' 2 

2 K ' 

2, 1 — k 



with ft = 0, 1, • • • , (N - 3)/2 and m = 1, 2,3, • • • . The prefactors 6 



N 
2 

- m 

N 
2 



;i 



(6.11b) 



(6.11c) 



1,2,3,4 to the 



hypergeometric function 3F2 are determined by the initial conditions ( 5. 10c ) and ( 6. lOd ), 



p(o,l) 

^0 



b[ L) - V2N 



p(o>2) 



2tt 



JV-2 \ 1 ' 



(6.12a) 



3 (o,2) 



(fc 

V^2 



^(2fe+l) 



3(0,2) 



fc+ i (JV-2fc-l)!! 

(Ar_l)l!(iV_2)!' 



A< 



U-F2 



1 _ 1. _fc _I 

2 5 5 2 

1 _ ^ 

2 ' 2 



(6.12b) 
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with k = 1,2, ■ •• ,(N- l)/2, and 



d(o,3) 



(fc + ^i)! 1 



& (2fc+2) _ g(°' 4 ). 



N 



h^2 



1 - h -k 1 

3 _1U _ JV 
2 ' 2 



;M) = ^t (^-2* -3)H 



(iV-3)!!(^2)!' 



(6.12c) 



with A; = 0, 1, • • • ,(iV-3)/2. 

In order to match this solution when 1 < m <C l/u> to the solution when N <C m, we need 
the leading and subleading terms of the large m behavior of (% ; and 6, n . They are 



with Cq "* some number that depends on TV, 



+ A k ~ 2 {N + 2k-l)(N - 1) 
with jfe = 1,2,- •• ,(iV- l)/2, and 

/ /o(o,4) \ 
h (2fc+2) _ 4 fc+i [ g(o,3) _^ _^fe \ m ^±i + fe 



>(o,D 



-+fc-l 



2fc + 2 



+ 4 fc_1 (JV + 2fc + 1)(JV - 3) ( e^ u ' o; + 



p(°,4) 
(o,3) . Ifc 



t+fc-1 



2fc + 2 / 

with = 0, 1, ■ • ■ , (N — 3)/2, to leading and subleading order in m. 



(6.13a) 



(6.13b) 



(6.13c) 



6.1.3. Matching solutions when N -C rn <C 1 

Having solved the recursion relations ([OJ) in the two limits N <C m and 1 <C m <C 
we are going to match the m = x 2 /(Alo) dependences of the solutions (5.9) and (5.13) in the 
overlapping region N <C m -C where both solutions are valid. We start with the coefficient 



bm . Matching Eq. ( 6.13a ) with Eq. ( |6.9a ) gives the two equations 



(0(0,1) _ 



c o,l £=i , /„(o) \ 



(6.14) 



for the two unknowns 6 



(i) 



and with the solutions 
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„(1) 



-26, 



(o) 



ft (D 26^ (V)! 

1 in 



(6.15) 



n/2/V. 



Matching Eq. ( p.l3b[ ) with Eq. ( p.9b| ) gives two equations (whereby we neglect In m and In Cj, 
compared to In uf) 



(o) 



-+k 



3(0,1) 



^0,1 



2(2fc)! 



In w = 4 te 2 



2fc 



(°,2) 



„(2fc+l) w 



-+fc-l 



3(0,1) 



(6.16a) 



4 k - 2 (N + 2k~ 1)(N -1) 



_£ i— P 

2k 



(0,2) 



o- 1 2(2fc-l)! 

for the two unknowns b^ k+1 ^ and c^ k+1 \ while matching Eq. ( |6.13c ) with Eq. (6.9c) gives two 
equations 



„(2fc+2) ^ 2 



2fc + 2 



J2fc+2) 



+fc-l 



— = 4*" 1 (N + 2k + 1) (N 3) I 6i°< 3) + A_ 



(°,4) 



2(2fc 



for the two unknowns f ) ^ 2fe+2 - ) and c^" ri "'. Solving Eqs. (5.16a) and (6.16b) for c l [ and b\ 
gives 

q$ = (6.17a) 

with i = 2, 3, • ■ ■ , AT, 



J2fe+2) 



,(') 



(6.16b) 



b (2k+l) = g(o,2) 



-2ft 



TV - 1 



+ ft ! + 



iV 



fc ! 3 ^2 



1 u u 2_ 

2 ft/, fv, g 

1 _u _ N 
•>.>"' 0. 



(6.17b) 



with k = 1, 2, • • • , (N - l)/2, and 



& (2fe+2) _ g(o,4) 



1 /A" + 1 



2fc + 2 



+ k ! + 



1 /7Y 



.1 _ z» -k I 
3 -h - M. 

2' * 2 



(6.17c) 



with ft = 0, 1, ■ • • , (iV — 3)/2. Observe that the coefficients (6.17a) of the modified Bessel 
functions K 2 ,K 4 ,- ■ ■ ,K N _ X are trivially independent of u, while the initial values (5.17b) 
and (5.17c) of the u> — recursion relation are independent of w in view of Eqs. (5.12b) 
and (6.12c), respectively. 



6.1.4. 10 dependence of the recursion relation ( 5.2) 

r — -1 r- (i) U) 

Equipped with the solutions ( p.lSD and ( p.l7| ) for the u dependence of b\ and c [ with i = 

(i) 

1, 2, • • ■ , N, we can derive the u dependence of the coefficients 6 r „ with m = 2, 3, • • • in the 
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limit 1 < m <C u 1 from Eq. (6.1 1 ). Insertion of Eq. (5.15 ) into the prefactors (5.12a) appearing 
in Eq. (5.11a) gives 



In 



(CPJ) to + ° (m-i)! 



(6.18a) 



with Cq ^ some n umber that depends on N. Insertion of Eq. ( 6.17b ) into the prefactor ( |6.12b ) 
appearing in Eq. ( 6.1 lb[ ) gives 



b (2k + i) _ e (o.2) (m + fc-l + f)! 



i — — A" i — 

O "■'1 "'1 n 



. 1 k 711/ 2 



fc-| + #)! 



(m - 1)! 



i — fc, 1 — fe, 1 — m 



(6.18b) 



with fc = 1, 2, • • • , (JV — l)/2. Insertion of Eq. (6.17c) into the prefactor (6.12c) appearing in 
Eq. (6.1 lc) gives 



h (2k+2) _ e (o,4) (m + fc-l + f)! 



(m-|)! 



i*2 



1 



1 — k — m 



p (o,4) 1 (m + fc-| + f)! 
fc 2fc + 2 (m-1)! 



j — k — m j 



(6.18c) 



with k = 0, 1, • • • , (N - 3)/2. 

In conclusion, the w dep e ndenc e of the solution to the recursion relation (6.2) follows from 
combining Eqs. ( |6.7| ), (6.8), ( 6.15 ), ( 6.17a ) when m ^S> N together with Eq. ( 6.18 ) when moj <C 
1. This gives 



-e, 



6 (1) = 



b (2k+l) = 



b (2k'+2) 



-to 2 K (2y/mu)) , 



1 (m-l-i+^1)! (0 ,2) (m-l + f)! 

ln ( C o 0) ") 1=0 ( 2J + 1 )(™-!-0' ( m -I)! = 



0, 



N < m, 



( |6.18b| ), m<w _ \ 

N < TO, 



0, 



( |6.18c| ), to«cj-\ 



TV < TO., 



TO <C W 



(6.19) 



withfc = 1,2, • • • , (7Y-l)/2an d k' = 0,1, ■ ■ ■ , (iV - 3) /2. Equation (gig) should be compared 
to Eq. (5.14). Observe that Eqs. (5.18b) and (5.18c) are independent of to. 



6.1.5. Leading energy dependence of the density of states 

We are ready to extract the leading behavior of the average density of states v{e), Eq. (3.25), 
for asymptotically small energies e. The density of states in the long wire limit Eq. (3.24) is given 
by 
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v{e) = jim is Tr-^e L (<p \ (B - B) |^ } R . (6.20) 

with the left and right ground states ( 4.10a ) and ( 4.10b| ). The normalizations (-k8) and ( |4.1l| ), the 
biorthogonal relations (p~9|), and the identities 

(B-B) |0) = JV|0), 



(B-B) 



,(2n+l) 
'R 

. (2n+2) 



(JV-4n-l) \m)£ n+1) + \m) L 



(2ra+l) 



(6.21) 



+ 2n+2) 



(B B) \m)^ = (N - An - 3) \m)^ |m) L 
deliver the expectation value 

oo (W-l)/2 2 oo (AT-3)/2 

l ^oI(b-s)i^o)h=^+E E ^ n (^ n+1) ) E E « 



,(2n+2) 



m— 1 n— 

(iV-l)/2 



?7i— 1 n— 



'(JV-l)/2 



AT- 



m— 1 n— 



E E E K^,i 

^ fc=0 
r -3)/2 

E K^, 



oo (JV-3)/2 

E E 



• & (2fe+l) 



m— 1 n— 



' ^(2fe+2) 



fc=0 



(6.22) 

The matrices M .jy x and M .jy 2 are given in Eq. (6T), while the normalization factors Nm]n 

(2) \ k & — -n 

and lim n are given by Eq. (4.8). In view of Eq. (5.19), it is necessary to break the sum over m 
in Eq. (5.22) into two parts separated by the integer m , with N <C m >C . We choose m 
to be some number independent of ui. With this choice and for sufficiently small u, we find that 
the sum over m in Eq. ( 6.22 ) is dominated by th e contributions from bm with m > m , i.e., by 
the contributions from the first line in Eq. ( |6.19 ), 



(JV-l)/2 



l (Vol (£ - 5) bo) R « / dm E ( 



(6.23) 



By inserting the solution for bm from Eq. ( |5.19 ) we obtain 



4fe (o) ) 2 7 

l (Vol (5 - 5) IVo)r « ; 2 \ ° (o : / dm§ i °\m)m N - 1 K 2 (2^) 

ln (Co Vw 



with the combinatorial factor 

(N-l)/2 

S^(m):= E ^!n(K]v,i 

71 = 

Using the substitution x 2 = 4mw, the integral transforms into 

4(e o) )V jV ' 



l (voIGb-b)|^ ) r 



In 2 (C* 



dx § 



ar 
4w 



2JV-1 



Kl (x) 



with .x := 2y/m uj. Let us now expand in small u> (i.e., large argument) 
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(6.24a) 



(6.24b) 



(6.25) 



where Sfy 1 ' is some number that only depends on N. In order to derive the asymptotic density 
of states, it is sufficient to keep the leading u> dependence in Eq. ( |5.26 ). This gives 



/ I PS d\ I \ ~ 4 ( e O) ) 2 Nuj 1 Q (o,l) 

(vol [B- B ) I<A)}r ~ 



In 2 (Cj 



(^f- — §V' / dx xK ( x ) 



(6.27) 



Combining Eq. ( J6.27 ) and ( 6.2C ), we find that the density of states, for asymptotically small 
energies, is given by 



v(e) oc lim Re ■ 



(o) 



e ln d £ 



(6.28) 



Recalling that e is measured in units of the disorder strength, we have recovered Eq. (1.1) when 
the number of channels N is odd. 



6.2. Density of states for an even number N of channels 



The calculation of the density of states for an even number N of channels follows along the 
lines of the N = 2 case with the caveat that for N > 2 the recursion relations (4.13) no longer 
decouple. This difference complicates the calculation considerably. To overcome this difficulty 

(n) 

we introduce a linear transformation of the coefficients a m that approximately decouples the 



recursion relations (4.13) in the two regimes m ^> N and m <C ui 1 . That is, we introduce 



N/2-1 

e +i) = e p*w]^ n+i) , 

N/2-1 



e +2) = e p^.d; 



k n (2n+2) 



N 

fc = 0,l,--- ,— -1, 



k = 0,l,--- -1, 



(6.29a) 



with m = 1, 2, 3, • • • and the transformation matrices 

» (2n)!(iV-l-2n)!(^) - 
L^iJn-l ^ -„)!(#-!)! 



-3F2 



-fc, fc + 1, — n 



1 2-iV 
2 ' 2 



r M , (2n+l)!(jV-2-2»)!(^)! 



F 9 



—A;, k + 1, — n 

3 2-jV 

2 ' 2 



1 



(6.29b) 



Here, the symbol 3F2 denotes a generalized hypergeometric function (see for example Ref. [ 93 ]) 

(n) (k) 
It is also possible to express the a) n s in terms of the b m s. Indeed, with the help of identi- 
ties (B.5a) and (B.5b), inverting the linear relation (5.29a), i.e., for any m = 1, 2, 3, ■ • • , 
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N/2-1 



a 



(2n+l) 



,(2ra+2) _ 



72-1 

J r -\ n 

/ j [ JVL o;Ar,2j "m 



(2fc+2) 



fc=0 



n = 0,l,.-.,y-l, 



N 

n = 0, 1,- • • , — - 1, 



(6.29c) 



can be done with the inverses 



M 



e;JV,l 



^,2 



{2n)\ 



[(N — l)!] 2 



fc (7V_i_2n)!(^-fc)!(^ + fc)! 

» = ^ I [(jV-2)!] 2 (2fc + l) 2 
fe ~ (JV-2-2n)!(^-*)!(^ + fc) 



3-F2 



;3-p2 



2-iV 
2 ' 

1 _ JV 2-N 

2 2 ' 2 



2-7V 



fc 2 , fc , 

-, _ JV 3-jV 
2 ' 2 



;i 



(6.29d) 



of the transformation matrices (5.29b). With these definitions, we are going to rewrite the re- 
cursion relations (4.13) in a form that decouples both in the limit c = and in the large m 
limit. 

To see this, we first make use of identities (B.3c) and ( ft. 3d ) to combine Eqs. (4.13) and (S. 29 ) 
into the recursion relations 

4mub£ k+1) = {-4m(2m - 1 + N) + [N - 2 - Ak(k + 1)]} b% k+1) 
+ 2m(2m + l)&S£ 1) 



+ (2m - 2 + N) (2m - 1 + N)b {2k _X 1] 
N/2-1 

-4c £ n[M c;7Vil ]^" +1 > 

71=0 



and 



4mc&< 2A;+2) = {-4m(2m - 3 + iV) + [3 AT - 6 - 4fc(fc + 1)]} 6^ fc+2) 

+ 2m(2m - l)b {2k ^ 2) 

+ (2m - 3 + AO (2m - 2 + N)b {2k \ 2) 
N/2-1 

-4c 2 n[M e;JV)2 ]^r +2) 

n=0 

with = 0, 1, ■ • ■ , (N/2) — 1 and m = 2, 3, • ■ • . Here, the initial conditions are 

4wb { i k+1} = + V2N(N + 1) - [3AT + 6 + Ak(k + 1)] b {2k+1) + 6b {2k+1) 



N/2-1 



4c n [W C ;N,i] n ai 



k (2n+l) 



(6.30a) 



(6.30b) 



(6.30c) 



with k = 0, 1, ■ • • , (N/2) - 1 and 

Acub[ 2k+2) = - V2N(N - 1) - [N + 2 + Ak(k + 1)] bf k+2) + 2b {2k+2) 



N/2-1 



4c J2 n [ M e;N, 2 ] n 4 



k (2n+2) 



(6.30d) 
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with k = 0, 1, 2, • ■ ■ , (N/2) — 1. We aim at a solution of Eq. (6.30) for asymptotically small 
frequency to. As in Sec. 5.2, our strategy is to solve Eq. (5.3C) in the two limits N m and 



m <C oj .In the former limit, m ^> N, we treat m as a continuous variable. Then, Eq. (5.3C) 



decouples, provided we assu me tha t a m decays rapidly for m J> N . In the other limit, m <C 
lo^ 1 , we can neglect w in Eq. ( 5.30 ), which again decouples Eq. ( 6.30 ). In this way, it is possible 
to find approximate solutions in the two regions N <C to and m <C lo^ 1 that are uniquely fixed 
up to some multiplicative coefficients and initial values, respectively. For asymptotically small 
u>, the overlap between these two regions N <C m <C is large. We can then match the two 
approximate solutions in the overlapping region. This gives a unique and approxi mate so lution 
for the ground state wave function, which in turn determines the density of states (3.25) in the 
long wire limit. 



6.2.1. Solution when m 3> N 

First, we treat the limit to 3> N. If we assu me that am decays rapidly for large to, we can 
drop the last line on the right-hand side of Eqs. ( |5.30a ) and ( |6.30b ), respectively. Consequently, 
in terms of the finite differences 



Abl*) := b. 



6 W 



A 2 &W := b 



m+2 



26. 



rn+1 



6 (0 

"TO i 



(6.31) 



with i = 1 , 



N, the recursion relation (6.30) reads 

(2fc+l)\ _ 



(m + l)4c(Ae +1) +e +1) ) 



2(2to + 2)(7V-2)A&0^ +1) 
(2to + 2)(2to + 3)A 2 ^+ 1 ) 
(AT - l)(AT-2)^ 2fc+1) 
[JV - 2 - 4As(fc + 1)] (A&^ 



(6.32a) 



j,(2fc+iA 



with A; = 0, 1, • • • ,(N/2)- 
(m + l)4w (A6, ( f'+ 2 ) 



1 and 

j. 6 (2fe+2) 



2(2TO + 2)(A^-2)A6^ fc+2) 

(2to + 2)(2to+ l)A 2 fe, (2fc+2) 
(iV-2)(^-3)fe, (2fe+2) 
[37V - 6 - 4fc(fc + 1)] (A&(f +2) 



W2fc+2) 

(6.32b) 

with k = , (AT/2) - 1. Observe that Eqs. ( |6.32a| ) and ( |5.32b| ) only differ from Eqs. ^A^) 



and ( |6.4b[ ) through the coefficients [N —2 — 4fc(fc + 1)] and [3N - 6 - 4fc(fc - 1)], respec- 
tively. In the limit m ^> N, we can neglect terms of order N compared to m and replace finite 



differences by derivatives. In place of Eq. (p.32|) and if we assume that b 
i = 1,2, • • • ,N is slowly varying, we get 

db (2k+l) 



(0 



& W (m) with 



4TOw& (2fc+1) = - Am(N _ 2)- 



dm 



4to 



,d 2 b( 2fc+1 ) 
dm 2 



+ [(N — l) 2 — (2fc + l) 2 ] fo( 2fe+1 ) 



and 

4TO^&< 2fc + 2 ) 



-4m(iV - 2) 



rffr( 2fc + 2 ) 
dm 



4to 2 



|d 2^(2fc+2) 

dm 2 



[(N 



(2k 



(6.33a) 

1)21 & (2fe+2) 

(6.33b) 
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with k = 0, 1, ■ • • , (N/2) - 1. By use of the substitution 



whereby 



4o> 



Aid T7l 



» N, 



we find that the solutions to Eq. (6.33) are given by the linear combinations 



2k+r 



(x) 



(2k+l) 



and 



b (2k+2 \x) = cf[ 



(2k+2) fX 



N-l 



N-l 



R- fr) 4- r (2fe+2) 



(f) 



AT-1 



L 2k+ 



(6.34a) 



(6.34b) 



(6.35a) 



(6.35b) 



with k = 0, 1, • ■ ■ , (N/2) — 1, of modified Bessel functions K 2k+1 and I 2 k+i- We shall demand 
that &m with i = 1, • • • , TV decay to zero as ?n — > oo, i.e., we must set 



$ = 



(6.36) 



with i = 1, • ■ ■ , N in Eq. ( p5| ) . The remaining iV coefficients c£\ with i = 1 , ■ • ■ , N of the 
modified Bessel functions ii^, j£T 3 , i^ 5 , • • • , K N _ 1 are fixed by matching solutions (5.35) to the 
solutions in the m <C uj^ 1 region. Thereto, we need to extract the terms that are of order x N+2k 
and x N+2k ~ 2 from the expansion of Eq. ( 6.35 ) when x small (see chapter 9.6 in Ref. [ 92 1). This 
gives 



b^ 2k+1 \x)^^ 2k+1) 



- c, 



(2fc+l) 



(x/2) 



N+2k 



and 



b^(x)~c?* +2) 



„(2fc+2) 



2(2k+ 1)! 

(mu>)~ +k 
2(2k + 1)! 

(a;/2) w + 2fc 



7 + In ■ 



2(2fc + 1)! 
2(2fc + 1) 



In (Cfmw 



7 + In ■ 



*(2fc + 2) 
(mw)^ +t_1 



(x/2) 



7V+2fc-2 



2(2fc)! 



- In (c^mu/) 



2(2fc)! 

*(2/« + 2) 
(mw)^ + 



(6.37a) 



(x/2) 



N+2k-2 



2(2k)\ 



k-1 



2(2k)\ 



(6.37b) 



with k = 0, 1, ■ • ■ , (N/2) — 1 when x <C 1, and where 



C 



(e) 



exp (7- *(2fc + 2)) 



(6.37c) 



6.2.2. Solution when 1 < m <C 

S econd , we treat the lim it 1 < m <C in which case we can neglect the ui terms in 

Eq. ( 6.30| ). In doing so, Eq. ( 6.30| ) becomes 



= {-4m(2m - 1 + N) + [IV - 2 - 4fc(fc + 1)]} 6^ fe+1) + 2m(2m + 1)6^ X) 



+ (2m - 2 + /V)(2m - 1 + iV^m-t^ 



(6.38a) 
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and 



= {-4m(2m - 3 + N) + [3N - 6 - 4fc(fc + 1)]} b<£ k+2) + 2to(2to - 1)6 



(2fc+2) 
'm+1 



+ (2m - 3 + JV)(2m - 2 + 7V)6^+ 2) 



(6.38b) 



with k = 0, 1, • • • , (N/2) — 1 and m = 2, 3, • • • , together with the initial conditions 

= +V2N(N + 1) - [37V + 6 + 4fc(fc + 1)] fe+1) + 6b ( 2 2k+1) (6.38c) 

and 

= -V2N(N — 1) — [N + 2 + Ak(k + 1)] b (2k+2) + 2b ( 2 2k+2) . (6.38d) 
The solution to the recursion relation ( |6.38 ) can be expressed in terms of generalized hypergeo- 
metric functions 



,(2fc+l) _p(e,l) (m + fc-l + 

~ k (m-1)! 



^-3-^2 



— fc, i — fc, 1 — to 

2 ;i 

|,l-jfc-m-f 



(e,2) 



(m-i)! 



-3-F2 



(6.39a) 



and 



,( 2 fc+2) _ p (c,3) (m + fc-l + f)! 

°m — L fc 1 TV, 3^2 



(to - 1)! 



-fc,-i 



/c, 1 — m 



— j 1 771 ~2~ 



: 1 



(c . 4) (TO + fc ~ 
L /c 7 3Ti 3-^2 



(to- |) 



— , 77 k , — TIT.' 



3 3 



3 



k to 2 



(6.39b) 



with & = 0, 1, ■ • ■ , (iV/2) - 1 and to = 2, 3, ■ • • . The prefactors e[ c,i) i = 1,2,3 , 4 to the 



hypergeometric function 3F2 are determined by the initial conditions ( 5.38c ) and ( 5.38df ) through 



,(e,l) 



a(e,2) 



(*■ 

V^2 



6 (2 fe+ l)_ e (e,2)^ ( 



7V+ 1 



fc+ i (JV-2fc-2)!! 
(iV-2)!!(^i)!' 



and 



3 (e,3) 



V^2 



(N-2)\\ (£l=£ 



! 3-F2 



I _1 _ _ W 

2 ' 2 2 



(6.40a) 



AT- 1 



! 3-F2 



-k I — i- I 



3 1 



AT 



(6.40b) 

with fc = 0, 1, ■ • ■ , (N/2) — 1. It is interesting to note that the explicit k dependences on the right- 
hand sides of Eqs. ( |5.39a ) and (6.39b) are related to those on the right-hand sides of Eqs. ( |6.1 lb ) 
and (|6.11cj) by letting k — > k- (1/2) and k — > k + (1/2), respectively. 
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In order to match this solution when 1 < m <C l/u> to the solution when JV « m, we need 
the leading and subleading terms of the large m behavior of bm +1 ^ and &m fc+2 ' ) . For b m ^ and 



l(2) 



we have 



b\}, 



N 
N 



2)6, 



(e,l) 



ive 



(o,2) 



^ - (e^' 3) + m* + f ((iV - 2)<t 3 > + (TV - 4)e<°' 4 >) m*" 1 , 



(6.41a) 



(6.41b) 



to leading and subleading order in m, respectively. In general we have for the odd-numbered 
coefficients 



L(2fc+1) 



2fc + 1 



e^ 2) mf + fc 



+ 4 fc_ 3(/V + 2fc)(jV- 1) 



2fc + 1 



<t 2) M 



-fe-i 



(6.41c) 



and for the even-numbered coefficients 



& (2fc+2) 



g(e,3) 



3 (e,4) 



2fc + 1 



+ 4 fc ^5(iV + 2fc)(iV-3) I e^ 3) + 



2fc + 1 



,#+fc-i 



(6.41d) 



with fc = 1, • • • , (N/2) — 1, to leading and subleading order in m. We emphasize that these 
two formula are only valid for k > 1. Observe that the large m subleading behavior of b m ^ and 
bm differ from those of bm k+1 ^ and bm k+2S> with k > 1, re spect ively. We note again that the 
explicit k dependences on th e right- hand sides o f Eqs. ( 5.41c ) and ( |6.41d ) are related to those on 
the right-hand sides of Eqs. ( |5.13b| ) and ( |5.13c[ ) by letting k -> k - (1/2) and k -> fc + (1/2), 
respectively. 



6.2.3. Matching solutions when N -C m <C cj -1 

Having solved the recursion relations ( 3.3C ) in the two limits TV -c m and 1 -c m -C 
we are going to match the m = .t 2 / (4o>) dependences of the solutions ( |6.37| ) and ( 6.41 ) in the 
overlapping region N <C m -C w ~ 1 , where both solutions are valid. We start with the coefficients 



bm and bm . Matching Eq. (5.37a) when k = with Eq. ( 6.41a| ) gives two equations in the two 



unknowns b^ and (whereby we neglect mm and lnC^ compared to In a/) 



I 1 ) 



■W 2 



iV 



(tv - 2)e^ e,1) + ive 



3 (e,2) 



(6.42a) 



,(2) 



Matching Eq. (6.37b) when k 

„( 2 ) 



with Eq. ( |5.41b[ ) gives two equations in the two unknowns 



and c, 



e,l> 



c (2) 



1 



: w 2 In w = e 



(e,3) 



3 (e,4) 



J 2 ) 



1 



iV _i 
1 J- 



iV r 



(A^-2)e{, e ' 3) + (^ 



(6.42b) 
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Solving for b±~\ b^\ c^\, and cfl, we determine their u> dependences, 



..(21 



and 



=V2(N + 1) 



bf ] = V2(N - 1) 



+4e< e - 2 W-f 

8-N(N -2)u>hLu>' 

-4e< M W-f 

8- N(N -2)wlnw' 

C^ 2) (f)! [8-N 2 ujlnu;] 

8- N(N -2)wlnw ' 
e ( eA) (f)\[8-N(N -4)u>]nu>] 



(6.43a) 



(6.43b) 



8 — N(N — 2)ujIiiuj 

respectively. For k > 1, matching Eq. ( |6.37e ) with Eq. ( |6.41c ) gives two equations for the two 



unknowns ^ 2fe+1 ) and <4 2 1 fe+1 ' ) (whereby we neglect In m and In c£ e ' compared to In cj) 



J2fc+1) W 



2(2fc + 1) 



In = 4' 



fe I L fc , p(e,2) 

2fc+l +Lfc 



(2fe+l)W 2+fc ] 



2(2fc)! 



A"-2(N + 2k)(N - 1) 



2fc + 1 



(6.44a) 



(e,2) 



while matching Eq. ( 6.37b ) with Eq. ( |5.41d| ) gives two equations for the two unknowns b{ 



and c 



(2fc+2) 



(2fc+2) 



e.l 



(2fc+2) ^ 2 



2(2fc + 1)! 



lnc=4 fc [ CP } ' Cfe 



(e,4) 



2k + 1 



c (^+2) j^Ll = 4 fc -f (7v + 2fc)(iv - 3) ( ei e - 3) + 1 



e (e,4) 



(6.44b) 



l e.l 



2(2fc) 



2k + 1 



with k = 1, 2, • • • , (AT/2) - 1. Solving for b^ and c^, we find 

(2fc+l) _„ (2fc+2) _ „ 

C i — u, c , — u, 



(6.45a) 



,(2fc+l) _«(e,2), 



(2Jfe+l) ( Z + k ) l + 



2 fN + 1 



k\l 3 F 2 



-k,-h-k,-l 



N 



1 -i - k - 

2 ' 2 2 



(6.45b) 



and 



b (2k+2) = g(e,4) 



-1 (N ,Y 1 /iV-1 



3 I _ I. _ N 
2 ' 2 2 



(6.45c) 

for k = 1,2, •■ ■ , (N/2) — 1. We note that the explicit fc dependences on the right-hand sides 
of Eq. ( |6.45b[ ) and Eq. ( |6.45c| ) are related to those on the right-hand sides of Eq. ( |6.17b[ ) and 
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Eq. (5.17c) by letting k — > k — (1/2) and k — > fc + (1/2), respectively. Observe that the coeffi- 

3, • ■ ■ , ^tv_i are trivially independent of 



cients ( 6.45a| ) of the modified Bessel functions K x , K 
, while the initial values ( 3.45b|) and ( |5.45c ) of the 
in view of Eqs. (|6.40a) and (|5.40b), respectively. 



recursion relation are independent of 



6.2.4. uj dependence of the recursion relation ( 5.3C) 



Equipped with the solutions (|6.43[) and (|6.45[) for the to dependence of b± and with 



(i) 

1, 2, • • • , TV, we can derive the u depe ndence of the coe fficien ts 6 r „ with m = 2, 3, 



in the limit 1 < to <C to 1 from Eq. ( J6.39 ). Insertion of Eq.( 5.43b ) into the prefactors ( 5.4C ) 
appearing in Eq. (6.39) when k = gives 



,(l)_p(e,2) 



ft (2) _p(e,4)_ 



! + £)! 



i + -V 

2r (4-iVwlnw) . 



4(m- 1)! 



(m-f)> 



e (^)(^_i±f)! (4+7Vcjlncj) 



4(m- 1)! 



(6.46a) 
(6.46b) 



where we have disregarded terms of order to and higher. Insertion of Eq. (5.45b) into the pref 



actor (5.40a) appearing in Eq. (5.39a) gives 



h (2k+i) _ p (e,2) (m + fc-^ + f)! 

+ *- 



1 7^* 7T") 

2 ' Q 1/6 2 



D(e,2) 



(2fc + 1)- 



(m-1)! 



-3-P2 



-fc, ^ — fc, 1 — m 



1 , 1 — fc — m 



2 



(6.46c) 



with fc = 1, 2, • ■ • , (N/2) — 1, Insertion of Eq. (6.45c) into the prefactor (5.40b) appearing in 
Eq. (6.39b) gives 



,(2fc+2) _p(e,4) (m + k- | + f) ! 



IV 



e (e,4) 



1 



-3F2 



k - 1 



2k + 1 



(to - 1)! 



-k i - 

' 2 

3 3 - i 
_ 2 ' 2 

2 / ' p 

3-^2 



A- 2 

' 2 



iV 

TO 2 



;i 



fc,i 



AT 

TO - - 



(6.46d) 



with fc = 1, 2, • ■ • , (N/2) — 1. We emphasize again that the explicit k dependences on the 
right-hand sides of Eq. (|6.46c ) and Eq. ( 5.46d[ ) are related to those on the right-hand sides of 
Eq. ( [6.18b| ) and Eq. ( |6.18c| ) by letting k^k - (1/2) and k ->■ k + (1/2), re spect ively. 

In conclusion, the to dependence of the solution to the recursion relation (6.30) follows from 
combining Eqs. (5.35), (6.36), (5.43), (5.45a) when m 3> N together with Eq. (6.46) when 
mto <C 1. This gives 
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b {2) = 



b (2k+l) = 



b (2k+2) 



■ -10; 



(e,2) 



-0 B-N(N -2)u,\nu> m ^ K l (2\M^) , 
g(e,2) ( (r^j+f)! _ (m-l + ff)! 



(e,4) 
8 



ivr7^r^ — i — 2 K-t (2Jmuj) 

— N{N — 2)ui\nuj 1 V V / 



-46, 



(e,4) ^ (m-| + f)! _ (m-l + f)! 



AT < to, 

TO <C 

iV < TO, 



-a i 



4(m-l)! 



(4 + iVw lnw) 



TO <C W 



-1 



(6.47) 



N <C TO, 



fl6.46cj ), to<w _1 , 

0, iV < TO, 

( |6.46^ ), to<cj- j , 



with k = 1 , 2, • ■ ■ , (iV/ 2) - 1. Equation ( |6.47[ ) should be compared to Eq. ( p.28[ ). Observe that 
Eqs. (6.46c) and ( |6.46d ) are independent of u>, 

6.2.5. Leading energy dependence of the density of states 



We are ready to extract the leading behavior of the average density of states z/( g), Eq . ( 3.25 ), 
for asymptotically small energies e. The density of states in the long wire limit Eq. ( 3.24 ) is given 
by 

u{e)= Urn. n- 1 Re^(<p \(B-B)\ ( p ) R , (6.48) 

u— * — ze 

with the left and right ground states (4.10a) and ( 4.10b[ ), respectively. The normalizations (4.8) 
and (4.1 1 ), the biorthogonal relations (4.9), and the identities 



(B - B) |0) = N |0) 
(B - B) \m\ 
(B-B) 



(2n+l) 
R 

(2n+2) _ 



/AT A 1 \ I \(2ri+l) , > (2n+l) 

[N — 4n — 1) |to)^ + |to)l 



(6.49) 



'R 



, q\ i \(2n+2) I i (2n+2) 

= (TV — 4n — 3) |to)^ ' — |to)l ; 



deliver the expectation value 



N/2-1 



N/2-1 



+ 2 



m—1 n— 
2V/2-1 



m—1 n— 



'N/2-1 



■■N ■ 



m—1 n— 



E E E K^,x 

\ fe=0 

72-1 

E 



00 JV/2-1 

E E >e 

m—1 n— 



(2fc+l) 



(2fc+2) 



(6.50) 



The matrices Mg.^ 1 and M c .jy 2 and their inverses are given in Eq. (6.29) while the nor- 



malization factors 'Nmln and Nm' n are given by Eq. (|4.8[). In order to evaluate Eq. (6.50) it is 
necessary to break the sum over to into two parts separated by the integer 1 <C to ~ uj~ 1 . This 
gives 

L {<p \(B-B)\<p ) K = N + S 1 + S i (6.51a) 



K 2 ) 
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with 



m N/2-1 m N/2-1 

Si :=£ £ *SiU i^ n+1) ) -EE 



(2n+2) 



rn—1 n—0 

oo N/2-1 



m— m n + l n—0 



m— 1 n—0 
2 



S * : = E E ^(KkJ/^) -3^n( 



M" 1 



N, 2 



'^ 2 



(6.51b) 
(6.51c) 



Here, we made use of the fact that in the large m limit the coefficients bm and bm all 
vanish except for bm and b$ . We are going to compute the sum S 1 for m < m Q and the sum 
S 2 for m > m separately. 

To compute the sum m < m , 5 l5 we first need to derive the leading lo dependence of the 
coefficients am™ +1 ' ) and flm" +2 '. Inserting Eq. (6.46a) and (5.46c) into Eq. (5.29c) yields 



,(2n+l) 



(-l)"V2(f 



(m + n)(m- 1)! (f -ra- l)!n! 



[1 + (m + n)w lnw] 



(6.52) 



where n = 0, 1, • ■ • , (N/2) - 1. Inserting Eq. ( |5.46b| ) and ( |6.46d| ) into Eq. ( |5.29c[ ) gives 

^ ^ (-l)V2(f +m - 1)! r , 

«, ( '" +2) = - : j J 1, 771V ^TTTl 1 + (m + n w In w , (6.53) 

(m + nj(m — 1)! — n — lj!n! 

where n = 0, 1, • ■ ■ , (N/2) — 1. With this, the normalization factors ( |4.8| ), and the choice m = 
uj" 1 , we find 

l/u) N/2-1 



*i = E E fe(<e +1) ) 



rn— 1 n—0 



(2n+2) 



2 ^ (| +m-l)! ^ (l + (m + nMn^) 
(m-1)! ^ (m + n)(f -n-1 

m— 1 v y n—0 v ; \ 2 



In! 



• m - fJ ! It -™- i) 1 \ n 2 



' N 



AT 



I)' 



(6.54) 



which is to be compared with the telescopic expansion (5.32). To compute the leading behavior 

(6.55) 



of the density of states it is sufficient to retain only the lowest order in lo. This gives 

Si ~ —N — const x to In 2 lo, 
to leading and subleading orders in lo. 

For the sum over m > m , S 2 , we obtain by inserting the definition of "M.~ N 1 , 2 , Nm] n , 
and Jvfm^n, as well as the solutions for bm and bm from Eq. ( |6.47 ) 

(e), , {AN) 2 ujm N - 1 K^(2y^U) 



m£? 0+ i [8-N(N-2)ojlnu,] 2 
with the combinatorial factor 

N/2-1 



2) 



iV,l 



JN, m,ri I ^0 



4) 



N.2 



(6.56a) 



(6.56b) 
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After converting the sum in S 2 into an integral, we find 

S 2 = 



m„+l 



W ^ [8- N(N -2)lo\tilo} 2 



(6.57) 



Using the substitution x 2 = 4muo the integral transforms into 



(4A0 2 w 



2, ,1-JV 



N(N -2)ojlnoj} 



OO 

(£) (!) 



2W-1 



(6.58) 



with xq = 2y/ (m Q + 1)lo ~ 0(1), since we chose m = lo 1 . Let us now expand S]y in small 
lo (i.e., large argument). This gives 



8«f" x 



2w 



7V-1 







2w 



N 



(6.59a) 



where the prefactor to the lo n 1 term 



(f)' 2 



AT/2-1 

E 



[(TV-^ir^Frf-f^n 



(7Y-2n- l)!(2n)! 
[ (jV - 2 )!]'( aFl [_^_ n ,|_^l]) s 
(JV-2n-2)!(2n+l)! 



(6.59b) 



= 0, 



is in fact vanishing. Here, we have introduced the Gauss hypergeometric function 2-F1. Hence, 
we conclude that the leading term of S 2 is (at most) of order w, and can therefore be neglected 
compared to S 1 , whose subleading w dependence is of order uilnu). The above reasoning also 
justifies the neglect of the m > m sum in Eq. (5.31 ). 

Combining Eq. (5.51) and (6.55), we find that the density of states, for asymptotically small 
energies, is given by 



v{e)= lim tt- 1 Rc(N + S 1 + S 2 ) 

u — > — is 

oc lim 7r _1 Rc (— lo In 2 lo) 

lj — > — is 

oc e |lne| . 



(6.60) 



Recalling that e is measured in units of the disorder strength, we have recovered Eq. (1.1) when 
the number of channels N is even. 



7. Discussion 



The asymptotic limit < e <C 1 for the density of states ( |l.l[ ) of a disordered quasi-one- 
dimensional quantum wire in the chiral-unitary symmetry class was first derived in Refs. [|82|] 
and [ p3| ] using the DMPK approach. This approach is based on finding the stationary solution 
of a forced diffusive process obeyed by the Lyapunov exponents of a transfer matrix, whereby 
the length L of the quasi-one-dimensional quantum wire plays the role of time and the forcing 
term is proportional to the frequency lo related to the energy e by the analytical continuation 
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uj — > —is. The DMPK approach is geometric, for the Lyapunov exponents are none but the 
radial coordinates of a smooth manifold determined by the symmetry class to which the disorder 
belongs. This symmetry requirement, the quasi-one-dimensionality, and the implicit weakness 
of the disorder determines in a unique way the forced diffusive process in the thick quantum 
wire limit TV — > oo [ |82| , p7| . In this paper, we have derived the density of states using the 
alternative superspin approach. 

The superspin approach relies on interpreting the disorder-average retarded Green function 
as a "quantum thermal average", whereby the role of temperature is played by the length L 
of the disordered quantum wire and, for a disordered q uasi- one-dimensional quantum wire, the 
quantum partition functio n is g iven by Eqs. (3JS) and ( 3.10 ). The quantum evolution, although 
unitary accor ding to Eq. ( 3.22 ), is here generated by the supersymmetric and non-Hermitian 
Hamiltonian (3.8). In the chiral-unitary symmetry class, the supersymmetric and non- Hermi tian 
Hamiltonian ( |3.8| ) simplifies to the supersymmetric and non-Hermitian Hamiltonian ( 3.12 ) so 
that it can be interpreted as a gl (1 1 1) quantum superspin Hamiltonian, for it reduces to a bilinear 
form in the generators of the superalgebra gl(2\2) and, furthermore, it commutes with a gl{\\l) 
sub-superalgebra of gl(2\2). In the thermodynamic limit L — > oo, the quantum statistical a verag e 
is solely determined by the non-degenerate right and left ground states of Hamiltonian ( |3.12| ), 
whose existence is guaranteed by supersymmetry. Thus, the computation of the density of states 
in the thermodynamic limit amounts to the construction o f the right and left ground states ( 4.1C ) 
for the supersymmetric and non-Hermitian Hamiltonian ( 3.12 ). 

Solving for these right and left ground states is difficult because the right and left Hilbert 
spaces on which the quantum superspin Hamiltonian is defined are infinite dimensional. For 
comparison, the irreducible Hilbert space of a single S't/(2)-spin is (2s + 1) -dimensional. For 
any given number TV of channels in the disordered quasi-one-dimensional quantum wire, the right 
and left ground states are constructed from two limiting solutions to the recursion relation (4.13) 

(n) 

for their expansion coefficients a m . Here, n = 1, • ■ ■ , TV is the channel index and m = 1, 2 , ■ • • 
is the basis index in the right Hilbert space. More precisely, the recursion relation ( 4.13 ) was 
solved independently in the large m 3> TV and in the small mui <C 1 limits, respectively, with 
the help of a change of basis in the right Hilbert space, i.e., by trading the expansion coefficients 
dm for the expansion coefficients bm > k = 1, ■ ■ ■ , TV. In the large m ^ TV limit, the recursion 
relation (4.13) reduces to TV modified Bessel equations whose leading solution for small u> is 
governed by (i) the modified Bessel function K when TV is odd or (ii) the modified Bessel 
function K 1 when TV is even. Evidently, this parity effect also characterizes the solution to the 
recursion relation in the limit mu> <C 1 provided JV<m< after both solutions have been 
matched in their overlapping regime TV <C m -C l/u> of validity. In turn, the same parity effect 
has a counterpart for the density of states for asymptotically small values of e, since, upon the 

(1) (k) 

analytical continuation u> — > — ie, (i) the coefficients b) n dominate over the coefficients b m with 
k = 2, • • ■ , TV in the density of states when TV is odd, while (ii) the coefficients bm and bm 
dominate over the coefficients bm with fc = 3, • • • , TV in the density of states when TV is even. 

In the superspin approach, the expansion coefficients b^ when TV is odd and bm and bm 
when TV is even play the role, in the DMPK approach, of the smallest Lyapunov exponent when 
TV is odd and the first two smallest Lyapunov exponents when TV is even, respectively. The only 
effect from the number TV of channels is to determine the size of the asymptotic regime e<1 for 
which Eq. ( |1.1[ ) holds through the choice g 2 oc v F /(N 2 £), where v F is the Fermi velocity in the 
clean limit and £ is the mean free path in the Born approximation, for the variance in Eq. (2.7). 
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Appendix A. Symmetries of the Hamiltonian 



In this appendix we will show that H 2 [as defined in Eq. ( 3.12 )] commutes with the eight 
generators of the Lie superalgebra gl (1 1 l)®gl (1| 1), while commutes with the four generators 
of the diagonal sub-algebra gl(l\l) C <?i(l|l) © <?Z(l|l).It thus follows that H is invariant under 
all transformations of (S^Sr) induced by the Lie supergroup Gi(l|l) generated by the Lie 
superalgebra gl(l\l). We refer to Ref. [Q] for an introduction on supermanifolds, to Ref. | ^ ] 
for an introduction on Lie superalgebra, and to Ref. [Q] for a dictionary of Lie superalgebras. 

The first step in the proof consists in rewriting H 2 as a quadratic form in the generators of the 
Lie superalgebra gl{2\2). To this end, define the 16 operators E ab with a, b = 1, • • • ,4 by: 
bosonic operators in the nonbar sector 

N N 

E xl = -B:=-Y J f}f i + \N, E^^-Q-^-^bth-lN, (A.l) 

i=l i=l 

fermionic operators in the nonbar sector 

N N 

E 13 ee +W + := + E bjfi , E 31 ee -W_ := - E fK (A.2) 

i=l i=l 

bosonic operators in the bar sector 

_ N _ N 

E 22 = ~B :=-EfifI + lN, E AA = -Q := + E h h \ ~ 5^. (A-.3) 

i=l i=l 

fermionic operators in the bar sector 

E 2i EE +W + := E kfl , E i2 EE-W_:=-t f z b\ , (A.4) 

i=l i=l 

bosonic operators in the mixed bar and nonbar sector 

N N 

E 12 = -A + := - E fJi, E 21 ee +A_ := + E flfl, 

i=l i=l 
N - N +-+ 

E 34 ee +D + := - E &A. E 43 ee +D_ := + E 

i=l i=l 

fermionic operators in the mixed bar and nonbar sector 

N N 



E 14 = -S + := + E fih E 41 ee +S_ := + E ffrl 

i=l i=l 

E 23 = -S + := - E b\f} , E 32 ee +S_ := + £ hi ■ 

i=i 



(A.5) 



(A.6) 



The 16 operators E ab with a, b = 1, • ■ ■ ,4 are either bosonic or fermionic. Hence, any E ab with 
a, b = 1, • ■ ■ ,4 can be assigned the degree (or grade) or 1 if they are bosonic or fermionic, 
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respectively, a mapping that we shall denote by dcg(E ab ). Our rational for the definition of the 
16 operators E ab with a, b = 1, ■ • ■ , 4 is that if we define 



0, if a = 1,2, 

1, if a = 3,4, 



then it follows that 



(A.7) 



(A.8) 



deg(£ ah ) = ([a] + [&])mod2. 
We can now define the 256 supercommutators 

[E ab , E cd ] := E ab E cd - (_!)**(*.») ^(E cd ) Ecd E ab , a, b, c, d = I, ■ ■ ■ , 4, (A.9) 



and verify that they define the Lie superalgebra of gl(2\2) in the standard basis [97|, i.e., 

^a6^cd]=^a d -(-l) ([Ql+[bl)([cl+[rfl) ^ ch , a, b, c,d = !,-••, 4. (A.10) 



The first step of the proof is then completed after verifying that H 2 , as given by Eqs. (3.12) and 
(3.8c), can be rewritten as 



Ho 



<h 



{A+f + 2BB - 2SS+ - 2W-W+ + {AJf 



2W-W+ + 2SS + - {D+f - 2QQ - (L>_) 2 



(A. 11) 



or alternatively 



H 



<-h 



( E 12) 2 + 2^11^22 + 2£; 41 £: 23 + 2^31^24 + ( E 2l) 2 



2EaoE 



"13 T 2£, 32^14 ~ (^34) 



2 ^33 -^44 



( E 4 3 ) 



Observe that 

H u =u(B-B + Q-Q) =u (E 22 - E u + E Ai - E 33 ) . 
In the second step, one verifies that H 2 commutes with the 8 operators 



(A. 12) 



(A.13) 



Si :— 7, {-E33 ~ E a + E 34 + ^43) > 



22 T" E 12 



E21): 



£■1 '■- 2 (+"^13 + E 24 - E U ~ E 2Z) > £ 1 : - 2 ( _ ^31 _ E 42 + E 41 + E 32) > 
S2 : = 2 (~ E 33 ~ E 44 - E 34 ~ E 43) > ^2 : = ^ (~ E 11 ~ ^22 _ E 12 ~ E 2l) > 



(A. 14) 



1 



1 



£ 2 : - 2 ( +E ™ + ^24 + E 14 + ^23) > £ 2 : - 2 ("^l ~~ ~~ ~~ ^2) ' 

that generate the Lie superalgebra gl(l\l) ® gl(l\l) with the only nonvanishing supercommuta- 
tors Mr 



[g 1 ,e±]=±£± [m 1; £±]=t£?, {£+ £r} = s 1 + w 1 , 
[s a ,e^]=±e± [m 2 ,£±] =t£ 2 ± , {4,£ 2 ~} = s 2 + m 2 
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At last, one verifies that the full Hamiltonian H = H u + H 2 commutes with the 4 operators 



N 



IB :=Jd + ?£ 2 = E 

i=l 
N 

Q:=Si + S 2 =E( fo K+^), 
i=i 

N 

i=l 
N 



(A. 16) 



that generate the Lie superalgebra gl(l\l) whose only nonvanishing supercommutators are 

[Q-S,W ± ] = ±2W±, {W + ,W_} = Q + S. (A.17) 

Thus, breaks the symmetry of H 2 generated by (1 1 1 ) © <^(1|1) down to the one generated 
by ff Z(l|l). 

Appendix B. Useful identities 

The generalized hypergeometric functions have the series expansion [B3| 



p r q 



Pi i ' ' ' , Pp 



E 



(B.la) 



where a lt ■ ■ ■ ,a p and ■ ■ ■ , (3 q are complex-valued parameters, z is complex-valued with 
magnitude less than unity, and 

r(a + n) 

(«)„ : ~ 



r(a) ' 

with r(z) the gamma function, denotes the Pochhammer symbol. 
We introduce the short-hand notations 

k . k . ■ tt< 



(B.lb) 



^i[/e,m,n] := 3 F 2 
F 2 [fc,TO, n] := 3 F 2 
F 3 [fc,m,ra] := 3 F 2 
F 4 [fc,m, n] := 3 F 2 



— fc, 2 + fc, — n 



-fc, fc + 1, — n 

-fc, fc + 1, — n 
I, -to, 



(B.2a) 



(B.2b) 



(B.2c) 



(B.2d) 
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for the family of hypergeometric functions labeled by the positive integers k, m, and n. They 
satisfy the identities 

= [4n(2n - 2m) - (2m - 4k 2 )] F 1 [k, m, n\ 

+2(2n + l)(m - n)F x [k, m, n + 1] + 2n(2m + 1 - 2n)F x [fc, to, n - 1], (B.3a) 
= {4n(2n - 2m) - [6m + 4 - 4(fc + l) 2 ]} F 2 [fc, m, nj 

+2(2n + 3)(m- n)F 2 [fc,m,n + 1] + 2n(2m + 3 - 2n)F 2 [k,m,n - 1], (B.3b) 
= {4n(2n - 2m - 1) - [2m - 4k(k + 1)]} F 3 [fc, m, n] 

+2(2n+ l)(m- n)F 3 [fc,m,n+ 1] + 2n(2m + 3 - 2n)F 3 [fc,m,n- 1], (B.3c) 
= {4n(2n + 1 - 2m) - [6m - Ak(k + 1)]} F 4 [fc, m, n] 

+ (2n + 3)(2m- 2n)F 4 [fc,m,n + 1] + 2n(2m + 1 - 2n)F i [k 1 m, n - 1], (B.3d) 

respectively. 

Furthermore, we introduce the notation 



F 3 1 [fc, m, n] 
-F^ 1 [fc, m, n] 



(m- fc)!(?7T, + fc + 1)! 
(2fc + l) 2 



3-F2 



— k — m — 1, fc — m, — n 



— fc — m — 1, fc — to, — n 



(to - fc)!(m + fc + 1)!' 
with the positive integers fc, m, and n. They satisfy the identities 

Y^F^[k,m,n]F,[k,m,n'] = {-l T ^ _A__^ 8n ^ 

k—0 

m 

fc=0 

for n, n' € {0, 1, 2, • ■ ■ , m} 



m! (2m + 1)! 

n|F#,m,n'] = (-l)"4 m " !( "! ~ ", )! ^ 

m!(2mj! 



(B.4a) 
(B.4b) 

(B.5a) 
(B.5b) 



References 

[I] P. W. Anderson, Phys. Rev. 109, 1492 (1958). 

[2] P. A. Lee and R. V. Ramakrishnan, Rev. Mod. Phys. 57, 287 (1985). 

[3] B. Kramer and A. MacKinnon, Rep. Prog. Phys. 56, 1469 (1993). 

[4] B. Kramer, T. Ohtsuki, and S. Kettemann, Phys. Rep. 417, 21 1 (2005). 

[5] F. Evers and A. D. Mirlin, Rev. Mod. Phys. 80, 1355 (2008). 

[6] J. Frohlich and T. Spencer, Commun. Math. Phys. 88, 151 (1983). 

[7] F. J. Dyson, J. Math. Phys. 6, 1 199 (1962). 

[8] J. J. M. Verbaarschot, Phys. Rev. Lett. 72, 2531 (1994). 

[9] M. R. Zirnbauer, J. Math. Phys. 37, 4986 (1996). 

[10] A. Altland and M. R. Zirnbauer, Phys. Rev. B 55, 1 142 (1997). 

[II] M. Caselle and W. Magnea, Phys. Rep. 394, 41 (2004). 

[12] P. Heinzner, A. Huckleberry, and M. R. Zirnbauer, Commun. Math. Phys. 257, 725 (2005). 



48 



[13] S. Helgason, Differential Geometry, Lie Groups, and Symmetric Spaces (Academic Press, San Diego, 1978). 

[14] Slightly abusing terminology, we also call a quantum critical point that separates two insulating phases (instead of 

separating a metallic from an insulating phase) a mobility edge in a problem of Anderson localization. 
[15] F. J. Wegner, Z. Physik B 35, 207 (1979). 

[16] S. Hikami, A. I. Larkin, and Y. Nagaoka, Prog. Theor. Phys. 63, 707 (1980). 
[17] L. Schafer and R. J. Wegner, Z. Physik B 38, 1 13 (1980). 

[18] K. B. Efetov, A. I. Larkin, and D. E. Khemlnitskii, Zh. Eksp. Teor. Fiz. 79, 1 120 (1980) [Sov. Phys. JETP 52, 568 
(1980)]. 

[19] A. Houghton, A. Jevicky, R. D. Kenway, and A. M. M. Pruisken, Phys. Rev. Lett. 45, 394 (1980). 
[20] S. Hikami, Phys. Rev. B 24, 2671 (1981). 
[21] K. B. Efetov, Adv. Phys. 32, 53 (1983). 

[22] H. Levine, S. B. Libby, and A. M. M. Pruisken, Phys. Rev. Lett. 51, 1915 (1983); A. M. M. Pruisken, Nucl. Phys. B 

235, 277 (1984). 
[23] R. Oppermann, Physica (Amsterdam) 167A, 301 (1990). 

[24] R. Gade, Nucl. Phys. B 398, 499 (1993); R. Gade and F. Wegner, ibid. 360, 213 (1991). 
[25] M. R. Zirnbauer, Phys. Rev. Lett. 69, 1584 (1992). 

[26] R. Bundschuh, C. Cassanello, D. Serban, and M. R. Zirnbauer, Nucl. Phys. B 532, 689 (1998); ibid Phys. Rev. B 
59, 4382 (1999). 

[27] T. Senthil, M. P. A. Fisher, L. Balents, and C. Nayak, Phys. Rev. Lett. 81, 4704 (1998); T. Senthil and M. P. A. 

fisher, Phys. Rev. B 60, 6893 (1999). 
[28] T. Senthil and M. P. A. Fisher, Phys. Rev. B 61, 9690 (2000). 
[29] M. Bocquet, D. Serban, and M. R. Zirnbauer, Nucl. Phys. B 578, 628 (2000). 
[30] N. Read and D. Green, Phys. Rev. B 61, 10267 (2000). 

[31] S. Guruswamy, A. LeClair, and A. W. W. Ludwig, Nucl. Phys. B 583, 475 (2000). 

[32] P. Fendley, Phys. Rev. B 63, 104429 (2001). 

[33] A. Altland and R. Merkt, Nucl. Phys. B 607, 5 1 1 (2001). 

[34] A. Lamacraft, B. D. Simons, and M. R. Zirnbauer, Phys. Rev. B 70, 075412 (2004). 
[35] S. Ryu, A. Furusaki, A. W. W. Ludwig, and C. Mudry, Nucl. Phys. B 780, 105 (2007). 
[36] P. M. Ostrovsky, I. V. Gornyi, and A. D. Mirlin, Phys. Rev. Lett. 98, 256801 (2007). 
[37] S. Ryu, C. Mudry, H. Obuse, and A. Furusaki, Phys. Rev. Lett. 99, 1 16601 (2007). 
[38] J. T. Chalker and P. D. Coddington, J. Phys. C 21, 2665 (1988). 

[39] D. K. K. Lee and J. T. Chalker, Phys. Rev. Lett. 72, 1510 (1994); D. K. K. Lee, J. T. Chalker, and D. Y. K. Ko, Phys. 

Rev. B 50, 5272(1994). 
[40] Y. B. Kim, A. Furusaki, and D. K. Lee, Phys. Rev. B 52, 16646 (1995). 
[41] J. T. Chalker and A. Dohmen, Phys. Rev. Lett. 75, 4496 (1995). 
[42] R. Klesse and R. Metzler, Europhys. Lett. 32, 229 (1995). 
[43] C.-M. Ho and J.T. Chalker, Phys. Rev. B 54, 8708 (1996). 
[44] S. Cho and M. P. A. Fisher, Phys. Rev. B 55, 1025 (1997). 
[45] V. Kagalovsky, B. Horovitz, and Y. Avishai, Phys. Rev. B 55, 7761 (1997). 
[46] R. Merkt, M. Janssen, and B. Huckestein, Phys. Rev. B 58, 4394 (1998). 
[47] V. Kagalovsky, B. Horovitz, Y. Avishai, and J.T. Chalker, Phys. Rev. Lett. 82, 3516 (1999). 

[48] J. T. Chalker, N. Read, V. Kagalovsky, B. Horovitz, Y. Avishai, and A. W. W. Ludwig, Phys. Rev. B 65, 012506 
(2001). 

[49] M. Bocquet and J. T. Chalker, Phys. Rev. B 67, 054204 (2003). 

[50] H. Obuse, A. Furusaki, S. Ryu, C. Mudry, Phys. Rev. B 76, 075301 (2007); ibid, Phys. Rev. B 78, 115301 (2008). 

[51] B. Huckestein, Rev. Mod. Phys. 67, 357 (1995). 

[52] N. Read, unpublished, 1991. 

[53] D.-H. Lee, Phys. Rev. B 50, 10788 (1994). 

[54] M. R. Zirnbauer, Ann. Phys. (Leipzig) 3, 513 (1994). 

[55] D.-H. Lee and Z. Wang, Philos. Mag. Lett. 73, 145 (1996). 

[56] J. Kondev and J.B. Marston, Nucl. Phys. B 497, 639 (1997). 

[57] M. R. Zirnbauer, J. Math. Phys. 38, 2007 (1997) and 2197 (1999). 

[58] J. B. Marston and S.-W. Tsai, Phys. Rev. Lett. 82, 4906 (1999). 

[59] I. Affleck, Nucl. Phys. B 265, 409 (1986). 

[60] I. A. Gruzberg, A. W. W. Ludwig, and N. Read, Phys. Rev. Lett. 82, 4524 (1999). 
[61] J. Cardy, Phys. Rev. Lett. 84, 3507 (2000). 



49 



[62] E. J. Beamond, J. Cardy, and J. T. Chalker, Phys. Rev. B 65, 214301 (2002). 
[63] A. D. Mirlin, F. Evers, and A. Mildenberger, J. Phys. A 36, 3255 (2003). 
[64] J. Cardy, Commun. Math. Phys. 258, 87 (2005). 
[65] F. J. Dyson, Phys. Rev. 92, 1331 (1953). 
[66] E. R. Smith, J. Phys. C 3, 1419 (1970). 

[67] G. Theodorou and M. H. Cohen, Phys. Rev. B 13, 4597 (1976). 
[68] L. Fleishman and D. C. Licciardello, J. Phys. C 10, L125 (1977). 

[69] A. A. Ovchinnikov and N. S. Erikhman, Zh. Eksp. Teor. Fiz. 73, 650 (1977) [Sov. Phys. JETP 46, 340 (1977)]. 
[70] V. L. Berezinskii, Zh. Eksp. Teor. Fiz. 65, 1251 (1973) [Sov. Phys. JETP 38, 620 (1974)]. 

[71] A. A. Gogolin and V. I. Melnikov, Zh. Eksp. Teor. Fiz. 73, 706 (1977) [Sov. Phys. JETP 46, 369 (1977)]; see also 

M. Steiner, M. Fabrizio, and A. O. Gogolin, Phys. Rev. B 57, 8290 (1998). 
[72] T. P. Eggarter and R Riedinger, Phys. Rev. B 18, 569 (1978). 
[73] A. D. Stone, and J. D. Joannopoulos, Phys. Rev. B 24, 3592 (1981). 
[74] E. Witten, Nucl. Phys. B 188, 513 (1981). 

[75] J. P. Bouchaud, A. Comtet, A. Georges, and P. Le Doussal, Ann. Phys. (N.Y.) 201, 285 (1990); A. Comtet, J. 

Desbois, and C. Monthus, Ann. Phys. (NY) 239, 312 (1995). 
[76] D. G. Shelton and A. M. Tsvelik, Phys. Rev. B 57, 14242 (1998). 
[77] L. Balents and M. P. A. Fisher, Phys. Rev. B 56, 12970 (1997). 
[78] M. Bocquet, Nucl. Phys. B 546, 621 (1999). 

[79] I. A. Gruzberg, N. Read, and A. W. W. Ludwig, Phys. Rev. B 63, 104422 (2001). 

[80] O. Motrunich, K. Damle, and D. A. Huse, Phys. Rev. B 65, 064206 (2002). 

[81] C. Mudry, S. Ryu, and A. Furusaki, Phys. Rev. B 67, 064202 (2003). 

[82] P. W. Brouwer, C. Mudry, and A. Furusaki, Phys. Rev. Lett. 84, 2913 (2000). 

[83] M. Titov, P. W. Brouwer, A. Furusaki, andC. Mudry, Phys. Rev. B 63, 235318 (2001). 

[84] O. N. Dorokhov, Pis'ma Zh. Eksp. Teor. Fiz. 36, 259 (1982) [JETP Letters 36, 318]; P. A. Mello, P. Pereyra, and N. 

Kumar, Ann. Phys. (NY) 181, 290 (1988). 
[85] C. Mudry, P. W. Brouwer, and A. Furusaki, Phys. Rev. B 59,13221 (1999). 
[86] A. Huffmann: J. Phys. A 23, 5733 (1990). 

[87] P. W. Brouwer, C. Mudry, and A. Furusaki, Nucl. Phys. B 565, 653 (2000); P. W. Brouwer, A. Furusaki, C. Mudry, 

and S. Ryu, BUTSURI 60, 935 (2005) (see cond-mat/051 1622 for english version). 
[88] P. W. Brouwer, A. Furusaki, I. A. Gruzberg, and C. Mudry, Phys. Rev. Lett. 85, 1064 (2000). 
[89] O. Motrunich, K. Damle, and D. A. Huse, Phys. Rev. B 63, 224204 (2001). 
[90] I. A. Gruzberg, N. Read, and S. Vishveshwara, Phys. Rev. B 71, 245124 (2005). 
[91] J. E. Bunder and R. H. McKenzie, Nucl. Phys. B 592, 445 (2001). 

[92] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, Dover, New- York (1965). 
[93] I. S. Gradshteyn and I. M. Ryzhik, Tables of Integrals, Series, and Products, 5th Ed., Academic Press, London 
(1994). 

[94] L. Frappat, A. Sciarrino, and P. Sorba, Dictionary on Lie Algebras and superalgebras, (Academic Press, London 
2000). 

[95] B. DeWitt, Supermanifolds, Second edition, (Cambridge University Press, New York 1992). 

[96] P. G. O. Freund, Supersymmetry, (Cambridge University Press, New York 1986). 

[97] Y.-Z. Zhang and M. D. Gould, J. Math. Phys. 46, 013505 (2005). 

[98] G. Gotz, T. Quella, and V. Schomerus, J. Algebra 312, 829 (2007). 

[99] L. Frappat, V. Hussin, and G. Rideau, J. Phys. A: Math Gen. 31, 4049 (1998). 



50 



